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Abstract: We reinterpret the JIMWLK/KLWMIJ evolution equation as the QCD Reggeon field 
theory (RFT). The basic ’’quantum Reggeon field” in this theory is the unitary matrix R which 
represents the single gluon scattering matrix. We discuss the peculiarities of the Hilbert space on 
which the RFT Hamiltonian acts. We develop a perturbative expansion in the RFT framework, 
and find several eigenstates of the zeroth order Hamiltonian. The zeroth order of this perturbation 
preserves the number of s - channel gluons. The eigenstates have a natural interpretation in terms 
of the t - channel exchanges. Studying the single s - channel gluon sector we find the eigenstates 
which include the reggeized gluon and five other colored Reggeons. In the two {s - channel) gluon 
sector we study only singlet color exchanges. We find five charge conjugation even states. The bound 
state of two reggeized gluons is the standard BFKL Pomeron. The intercepts of the other Pomerons 
in the large N limit are l + wp = 1 f- 2uj where 1 + w is the intercept of the BFKL Pomeron, but 
their coupling in perturbation theory is suppressed by at least 1/V^ relative to the double BFKL 
Pomeron exchange. For the [27, 27] Pomeron we find cu[ 27 , 27 ] = 2uj-\-0{l/N) > 2 uj. We also hnd three 
charge conjugation odd exchanges, one of which is the unit intercept Bartels-Lipatov-Vacca Odderon, 
while another one has an interecept greater than unity. We explain in what sense our calculation 
goes beyond the standard BFKL/BKP calculation. We make additional comments and discuss open 
questions in our approach. 
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1. Introduction 

The theoretical approach to high energy scattering in the pre-QCD era was established by V. Gribov 
and is known as Gribov‘s Reggeon Calculus Q. The idea was further developed and formulated 
as an effective Reggeon field theory (RFT) in Refs. ii- Since the advent of QCD more than 
three decades ago great effort has been made to reformulate these ideas in the framework of the 


- 1 - 





















QCD Lagrangian and thus to obtain RFT from first principles as a bona fide high energy limit of 
the theory of strong interactions. Despite significant progress this goal has not been achieved yet. 
Although various elements of RFT in QCD have been available for some time, a coherent formulation 
of RFT is still not at hand. 


The study of high energy limit in QCD began with the derivation of the Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) Pomeron Q, which was obtained as an infinite sum of perturbative diagrams com¬ 
puted in the leading logarithmic approximation (LLA). The key element of the BFKL theory is the 
concept of gluon reggeization, which is proven both in leading order (LO) |^, and next to leading 
order (NLO) 0. In particular, the BFKL Pomeron is understood as a ’’bound state” of two reggeized 
gluons in the t-channel. 


The BFKL Pomeron has an intercept greater than unity and leads to scattering amplitudes that 
grow as a power of energy, hence violating the s-channel unitarity. It was shortly realized that in 
order to restore unitarity, which undoubtedly has to be satished by the full QCD amplitudes, one 
needs to consider t-channel states with more than just two reggeized gluons. This was put forward 
by Bartels and is today known as generalized leading log approximation (GLLA). The important 
milestone result was a derivation of the BKP equation, which governs (perturbative) high energy 
evolution of an amplitude due to an exchange of an arbitrary but fixed number of reggeized gluons 
m. It was further realized that the s-channel unitarity cannot be achieved without introducing 
transitions between states with different number of reggeized gluons Hi- This is known as extended 
GLLA (EGLLA) |11, |l^]. The vigorous program of solving the BKP equations in the large Nc limit 


has been pursued in recent years by Lipatov and Korchemsky with collaborators [1R|,| 

The ideas put forward in have been under intense investigation during the last decade or so 


113, |1|, |I|, ^ n, 

in 0, [^, [^] and 2 


6 in 


25]. The transition vertex between 2 and 4 reggeized gluons was derived 
Oj. At present, this approach therefore gives elements of an effective 


theory in terms of f-channel gluon states and transition vertices. These elements have been put 
together in by Braun into an effective theory of BFKL Pomerons interacting via triple pomeron 
vertex of | pT| , . This model is meant to describe nucleus-nucleus collisions at high energies at LLA 

and large N^- In principle it contains Pomeron loops even though they were not included in the 
semiclassical analysis of |^|. 

A different approach to high energy scattering in QGD has been developed in the last 10 years 
or so based on the ideas of gluon saturation. The seminal paper of Gribov, Levin and Ryskin [^] 
put forward the idea that nonlinear effects in QCD evolution lead to saturation of gluonic density at 
high enough energy. This gluon saturation than should restore the s-channel unitarity which is badly 
violated by the BFKL Pomeron. The GLR equation - the nonlinear evolution equation for gluon 
density in the double logarithmic approximation was derived in |27|. The gluon saturation ideas have 
been further developed in a series of papers by Mueller]^], who also introduced the notion of QCD 
dipoles as a convenient basis for the discussion of high energy processes, and has related it to BFKL 
Pomeron. The dipole version of the triple Pomeron vertex was obtained in Ref. |30||. 


McLerran and Venugopalan recast the problem of saturation into that of studying nonlinear- 
ities of the classical Yang-Mills field theory directly in the path integral approach [^I| . This pro¬ 
vided an impetus for the formulation of the nonlinear QCD evolution approach to high energy 
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scattering[^|,|^],|^|,|^|. This line of research lead to two equivalent formulations of the nonlinear 
QCD evolution. One is formulated in terms of infinite hierarchy of evolution equations for mul- 
tiparton scattering amplitudes!321 - the so called Balitsky hierarchy; while the other one is based 
on a functional evolution equation for the probability distribution functional that determines the 
probabilities of various field configurations in the hadronic wave function - the so called JIMWLK 
equation 1 34],^^. The relation between the two formulations is roughly the same as between the 
infinite set of Dyson-Schwinger equations and the second quantized Hamiltonian in a quantum field 
theory. In this paper we find it more convenient to use the JIMWLK formulation, and will refer to 
the functional evolution kernel as the JIMWLK Hamiltonian. 


In the last year or so further progress has been made triggered by the realization that the 
JIMWLK equation does not include Pomeron loops | 


37, 


39| . It was realized that the JIMWLK 
Hamiltonian must be extended in such a way that makes it self dual under the Dense-Dilute Du¬ 
ality transformation [40, 41 j. Although the complete selfdual kernel has not yet been derived, the 
dual of the JIMWLK Hamiltonian (which we will refer to as KLWMIJ) was shown to describe the 
perturbative evolution of a dilute hadronic wave function [^]. The efforts to consistently include 
Pomeron loops into the evolution are ongoing and some progress in this respect has been made 


1 45, ^ 46, 47, 41, 48 


I 51, 52, 


We note that there also has been earlier work on the subject of the high energy evolution which 
is similar in spirit. In particular |53] considered an effective action for high energy QCD with Wilson 
lines as effective degrees of freedom. The very same degrees of freedom appear in more recent 
studies of Refs. [Q, 49, Mj. Yet another approach is due to Lipatov and collaborators [|^], who 

derived an effective action with both real and reggeized gluons as effective degrees of freedom. This 
action respects the unitarity of full QCD, but its complexity has so far precluded any progress in 
understanding its physical consequences. 


Obviously establishing the direct relation between the JIMWLK formalism, which is firmly based 
in perturbative QCD and the Gribov Reggeon field theory ideas should be useful in unifying the 
two approaches and bringing closer the two ’’wings” of the high energy QCD community. Some 
attempts in this direction have been made starting with the classic paper of Gribov, Levin and 
Ryskin |^, where the Pomeron fan diagrams were analyzed. Most of the attempts in recent years 
have concentrated on the mean field approximation to the dipole (large Nc) limit of high energy 


evolution 


|, which is technically a much simpler problem. These approaches converge 


which sums fan-type Pomeron diagrams |]^ in LLA 

5| . Beyond the 


to the Balitsky-Kovchegov (BK) equation 

and describes the nonlinear corrections to the BFKL Pomeronj^j and the Odderon 
large Nf. limit it is known that the leading perturbative order of the JIMWLK evolution reproduces 
the BKP hierarchy |34| and the perturbative Odderon solution [^]. 


Still the relations studied so far are somewhat indirect and do not provide a complete picture 
for translating between the JIMWLK and the Reggeon languages. The aim of the present paper is 
precisely to hll in this gap. We demonstrate that the JIMWLK/KLWMIJ Hamiltonian is nothing 
but the proper quantum field theoretical formulation of the Reggeon field theory. Our approach in 
the present paper will be based on the KLWMIJ form of the evolution, but as will become clear later 
on using JIMWLK in this context is entirely equivalent. We note, that we do not address here the 
problem of Pomeron loops in the JIMWLK framework, and therefore the Reggeon field theory we 
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deal with is not the complete one. However the mapping between the two languages is independent 
of the exact form of the Hamiltonian and should also encompass the complete and as yet unknown 
self dual generalization. 

This paper is structured as follows. In Section 2 we briefly review the general setup of calculating 
scattering amplitudes in the high energy eikonal approximation, the JIMWLK/KLWMIJ evolution 
kernel and discuss the structure of the Hilbert space on which the KLWMIJ Hamiltonian acts. 


In Section 3 as a warmup to the main part of the paper we discuss the dipole limit of the JIMWLK 
evolution[47, 411. We show that when the dipole model Hamiltonian is treated as a quantum field 
theory, one can set up a formal perturbative expansion around its quadratic part. The quadratic 
part preserves the dipole number, and so perturbative solution can be given in a sector with fixed 
number of dipoles. Solving the quadratic part itself one finds eigenfunctions which correspond to 
multi - BFKL Pomeron states, with the eigenvalues given by the multiples of the BFKL trajectory. 
The process of the diagonalization can be viewed (with some disclaimers) as finding such s - channel 
dipole impact factors which couple exclusively to the n-Pomeron states in the f-channel. 


Section 4 is the main part of this paper. Our aim here is to study the 2+1 dimensional quantum 
field theory defined by the full KLWMIJ Hamiltonian and to relate its perturbative spectrum with 
the results of the perturbative BFKL/BKP approach. This is of course much more complicated than 
for the dipole model, as the basic degrees of freedom now carry color indices. In fact the symmetry 
of this Hamiltonian is SU{N) (g) SU{N) ® Z2® Z2. One discrete Z2 factor corresponds to the charge 
conjugation symmetry of QCD (C-parity). The second Z 2 appears due to the symmetry between s 
and u channel scattering amplitudes in the eikonal approximation. It naturally corresponds with the 
Reggeon signature extensively discussed in the literature. The SU{N) (g) SU{N) symmetry also has 
a natural interpretation as an independent rotation of the incoming and outgoing color indices of the 
gluon scattering matrix. In the perturbative regime, where the gluon scattering matrix is close to 
unity, this symmetry is spontaneously broken down to diagonal SUv{N). The disordered regime on 
the other hand corresponds to the black body limit, since in physical terms it means that the color of 
the outgoing state after scattering is completely uncorrelated with the color of the incoming one. Like 
in the dipole model we develop perturbation theory around a ’’perturbative vacuum state”, that is 
the state in which the single gluon scattering matrix is close to unity. We develop an approximation 
which we will refer to as the partonic approximation to the KLWMIJ Hamiltonian, which preserves 
the number of s channel gluons throughout the evolution. We split the KLWMIJ Hamiltonian into the 
perturbative part and the interaction part. The perturbative part can be solved separately in every 
sector with fixed number of s - channel gluons. However its structure is much more complicated than 
in the dipole model, since it is not quadratic but rather quartic in the gluon creation and annihilation 
operators. 


We solve the one gluon sector of the theory and find states which directly correspond to reggeized 
gluon as well as other nonsinglet Reggeons. In our setup the color quantum number of the state 
’’exchanged in the t channel” is easily identified as the representation of the eigenstate with respect 
to the unbroken SUv{N) symmetry. We find the octet d-Reggeon which is degenerate with the 
reggeized gluon as well as Reggeons corresponding to other irreducible representations which can be 
obtained from multiplying two adjoints 8® 8 = 1 + 8^ + 8^ + 10 + 10 + 27 + TZj The reggeized 

^We follow in this paper the notations and nomenclature of j^. Thus the different representations of the SU{N) 
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gluon wave function and trajectory are the same as in the BFKL calculation. The non octet Reggeons 
have all the same spatial wave function and their trajectories are simple factors (proportional to the 
second Casimir of the representation) multiplying the trajectory of the reggeized gluon. Notably the 
singlet Reggeon exchange leads to energy independent contribution to cross section. 


We then go on to solve the two particle sector, but limit our solution to singlet exchanges in the 
t-channel. The structure we find can be easily understood in term of ’’bound states” of the Reggeons 
we found in the single s - channel gluon sector. We find states which are even and odd under the 
charge conjugation and the Z 2 signature symmetry. The C-parity and signature even sector contains 
the bound state of two reggeized gluons - the standard BFKL Pomeron with the usual wave function 
and trajectory. The other C-parity and signature even states include the bound state of two octet 
d-Reggeons (which is degenerate with the BFKL Pomeron), as well as bound states of two 27 , two 
77.7 Reggeons and a (color)symmetric bound state of 10 and 10 Reggeons. All of those in the large 
Nc limit have an intercept l + wp = l + 2a; where 1 + a; is the intercept of the BFKL Pomeron. The 
1/Nc correction to the intercept of TT-y Reggeon bound state is negative. The 1/W correction to the 
intercept of the 27 Reggeon bound state is positive, making it the dominant exchange. We show that 
the [ 85 , 85 ], [27,27] and [ 77 . 7 , 77 . 7 ] Pomerons contain at least four gluons in the t-channel, while the 
[10 + 10,10 + 10] Pomeron contains at least six t - channel gluons, making its coupling subleading in 
the dilute target limit. 


In the (P-parity odd sector we find three odderon states (for a recent review on Odderon see 
ref. fS^ ). One of them is signature odd and corresponds to the bound state of the reggeized gluon 
and the d-Reggeon antisymmetric in color space and coordinate space separately. We show that 
the wave function of this state contains at least three gluons in the t - channel and it is therefore 
identified with the Bartels-Lipatov-Vacca Odderon (BLV) ]59j. Another Odderon state is signature 
even. Its interpretation is that of the antisymmetric bound state of the 10 and 10 Reggeons. Its 
wave function contains at least six t-channel gluons. Correspondingly its trajectory is twice that of 
the BLV Odderon. The intercept of both of these states is unity. The third state we find is quite 
peculiar and interesting. It has a negative signature and can be thought of as a bound state of the 
reggeized gluon and the d-Reggeon symmetric in color space and coordinate space separately. Its 
intercept is greater than one and equal to that of the BFKL Pomeron. In general this state has at 
least three t-channel gluons and thus its coupling is not suppressed relative to the BLV Odderon. We 
show however that it decouples from a quark dipole. 


In Section 5 we discuss the relation of our calculation to the BFKL/BKP setup. The BKP 
approximation arises as the limit of the partonic approximation to the KLWMIJ Hamiltonian when 
each gluon emitted in the process of the evolution is allowed to interact with the target only once. The 
partonic approximation is more general since it allows multiple scatterings of the emitted gluons. We 
also discuss the way in which the ’’composite states” of many Reggeons contribute to the unitarization 
of the scattering amplitudes and the limitations of this picture. 


We conclude in Section 6 with discussion of our results and open questions. Three Appendices 
contain details of our calculations as well as a summary of useful properties of projectors for SU(N) 
group. 

group are labeled by the dimensionality of their counterparts in SU{3), except for TZr which does not exist in the SU{3) 
case. 


- 5 - 




2. High energy scattering: the general setup 


The process we are interested in is the scattering of a highly energetic left moving projectile consisting 
of gluons on a hadronic target. We are working in the light cone gauge natural to the projectile wave 
function, A~ = 0. In this gauge the high energy scattering matrix of a single gluon projectile at 
transverse position x on the target is given by the eikonal factor ^ 


5(a:;0, X ) = P exp{i / dy T°'a^{x,y )}; 


S{x) = S{x] 0,1) 


( 2 . 1 ) 


where is the generator of the SU{N) group in the adjoint representation. 

The field ot is the large A'^ component created by the target color charge density. It obeys 
the classical equation of motion and is determined by the color charge density of the target pt{x) 
via@, g, ^ 

=g^ p^y,x~)T^ 5'(y;0,x")| (2.2) 


For a composite projectile which has some distribution of gluons in its wave function the eikonal 
factor can be written in the form analogous to S{x), see 


= / Dpp W^[pp\ exp 




d'^xpp{x,y )ap{x,y ) 


(2.3) 


with Xi - the transverse coordinate. The quantity pp{xi) is the color charge density in the projectile 
wave function at a given transverse position, while VF^[p] can be thought of as the weight functional 
which determines the probability density to find a given configuration of color charges in the projectile. 
For a single gluon p°'{xi) = T^5'^{xi — x°), and eq. (|2.3| ) reduces to eq. O- 

The total S'-matrix of the scattering process at a given rapidity Y is given by 


^(T) = j Dap WY^[aT{x,x )] S^_yjar(x,x )] 


(2.4) 


In Eq. (2.4) we have restored the rapidity variable and have chosen the frame where the target has 
rapidity Yq while the projectile carries the rest of the total rapidity Y. Lorentz invariance requires 
S to be independent of To- 

The high energy evolution is generally given by the following expression: 


d 

dY 


S= Da^ W^^[aT{x,x )] y 


ap, 


S ap 


Y^_YJapix,x )]. 


(2.5) 


Here y stands for a generic Hermitian kernel of high energy evolution, which can be viewed as acting 
either to the right or as Hermitian conjugated to the left: 


— 'pP- 
dY^ ^ 


ap, 


6ap 




ap, 


6ap 


W^[ap ]. 


( 2 . 6 ) 


^In our convention the variable x is rescaled to run from 0 to 1. 
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As was shown in |^] in order for the total 5-matrix to be Lorentz invariant and symmetric between 
the projectile and the target, the evolution kernel y must be self dual. That is it has to be invariant 
under the Dense-Dilute Duality transformation 


a“(x ,X 


5p°-{x, 


X 


^ ■ at 

- -— ^ -tp {x,x ) 

da^-lXjX ) 


(2.7) 


However if one considers the situation where the target is large and the projectile is small, the 
symmetry between the target and the projectile is irrelevant. In the limit when the color charge 
density of the target is parametrically large (/9“ = 0{l/as)) the kernel is given by the JIMWLK 
expression |^, ^] 


^JIMWLK _ f> 

X — -^x,y,z 


2 tr 


6Sl 




Sftr 


SyT^^ 

6Sl 



r 


r 


r 


r 


tr 

Lj51 j 

tr 

-^T'^Sy 

6Sl 

— tr 

SxT^^ 

1 J5lJ 

tr 

SyT^^ 

6Sl 

\ 


( 2 . 8 ) 




as {z-x)i{z-y)i 


^,y,z 27]-2 — xY {z — uY ' 

The functional derivatives in this expression do not act on in the kernel. 


Kx,y,zf = [ Kx,y,z f{x,y,z) (2.9) 

Jx,y,z 


Referring back to eq. (|2.3| ) we see that S^[a] is the functional Fourier transform of IT^[/5]. It 
thus follows that if evolves according to eq. (|2.6| )with the JIMWLK kernel, must evolve 

with the dual kernel, which we call KLWMIJ 




KLWMIJ 


^1 r 

5p_ 


W^[p]. 


( 2 . 10 ) 


with Q 


^KLWMIJ ^ J 2tr 


6Ri 


r“ Rx 


Rftr 


R„ T 




— tr 


[6R: 


T“ Rr 


tr 


6R 


T a D 

Ky 

I ^ 


— tr 


Rr r“ 


6Rl 


tr 


Ry r“ 


6rI 


( 2 . 11 ) 


where the ’’dual Wilson line” R is defined as 

R{zY^ = 


/■i _ S 

Rexp / dz 7 -ry 

Jo [z, z j 


ab 


( 2 . 12 ) 


Just like in Eq. (|2.8| ) the functional derivatives with respect to R do not act on R{z). This can be 
rewritten in the normal ordered form, that is with all derivatives acting only on external factors. 
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KLWMIJ 

X = -jXx,y,z ■ 


— 2tr 


6RI 


r“ R, 


Rftr 


R.T 


6rI 


+ tr 


6RI 


r“ R, 


tr 


6R. 


T a p 

Ky 


+ tr 


RxT^ 


6RL 


tr 


D n~'0- 
Uy I - T 

6Rl 


- Kx 


- 2 Rf tr 


r“ Rr 


SRl 


+ 2N tr 


SRl 


Rx 


(2.13) 


We will find it convenient to rewrite the KLWMIJ kernel in the following form 


KLWMIJ 

X = - A, 


x,y,z 


+ tr 


6R: 


T\Rx - R,) 


— 2tr 

6 


■T^{Rx - Rz) 


tr 


SR 


I SR] 


T%Ry - R, 


Rftr 


{Ry - 


sr\ 


+ tr 


{Rx - Rz)T° 


srI. 


tr 


{Ry - R,)T^ 


SRl 


-K, 


X^X,Z 


- 2 Rf tr 


T^{Rx - Rz)T^-^ 
SRli 


+ 2 A^tr 


ISRi 


{Rx - Rz 


(2.14) 


To arrive at this form we used explicitly the fact that J? is a unitary matrix in the adjoint repre¬ 
sentation. This form of the kernel is particularly convenient, since each term in eq.( |2.14 ) separately 
is ultraviolet finite. The complete KLWMIJ kernel is ultraviolet finite, that is apparent ultraviolet 
divergencies at the point x = y, z = x cancel between different terms in eq.( |2.13| ). In the form 
eq.( 2.14 ) this finiteness is explicit due to the factors R^ — Rx etc. 


Note that the KLWMIJ kernel possesses the global SUl{N)0SUfi,{N) symmetry. The symmetry 
transformations rotate the left and right indices of the matrix R independently. Furthermore there 
is a discrete Z 2 symmetry ^ This transformation interchanges the color states of the 

incoming and outgoing gluons. Thus this symmetry reflects the s — u crossing invariance of the QCD 
scattering amplitudes. The quantum number (charge) associated with the Z 2 symmetry goes under 
the name of signature. The fact that the KLWMIJ Hamiltonian preserves this symmetry is nothing 
else but the Gribov‘s signature conservation rule. Below we will see the signature conservation rule 
at work. The kernel is also charge conjugation invariant. We will discuss the Z 2 charge conjugation 
transformation in detail in Section 4. 


An important point which we have not mentioned so far, is that W[p\ can not be an arbitrary 
functional. This issue has been discussed in detail in [44, 48| and more recently in [61, 62|. The point 
is that the longitudinal coordinate x~ is introduced in our formalism as an ’’ordering” coordinate. To 
begin with the transverse charge density p°‘{x) depends only on transverse coordinates. However p°‘{x) 
are quantum operators and therefore do not commute between themselves. To emulate calculation of 
correlation functions of a product of such operators in the projectile wave function we have introduced 
in 1^1 an additional variable x~ whose values simply track the order of appearance of the operator 
factors of p°‘{x) in the product. The consistency requirement on W\p] is then such that classical 
correlation functions calculated with W as the weight functional must reproduce the relations between 
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the correlators of p°‘{x) which follow from the SU{N) algebra satisfied by p's. We have shown in 
that the functionals of the following form satisfy all the required relations 


W[p] = T.[R]5[p{x,x )] 


(2.15) 


with an arbitrary functional S which depends on the ’’dual Wilson line” R only. Eq.( ^.15 ) is a simple 
restriction on the Fourier transform of W. Inverting it we find 


J d/9exp|i J dx d^xp°'{x,x )a“(x,x )| W[p] = S[5] 


(2.16) 


with S defined as in eq.(2.1). Thus the Fourier transform of W must depend only on the Wilson 
line S and not on any other combination of the conjugate variables a. This of course has a very 
transparent physical interpretation. The functional Fourier transform of W is precisely the projectile 
averaged scattering matrix for scattering on the ’’external field” a. The restriction eq.( p.l5 ) simply 
means that this scattering matrix must be a function of scattering matrices of individual gluons 
constituting the projectile, and does not depend on any other property of the external field. The 
normalization of the functional W is determined by requiring that for a = 0 the scattering matrix is 
equal to unity. Thus 

F(i2 = l) = l, j dpW[p] = 1 (2.17) 

Further general properties of S follow from its identification as the projectile averaged scattering 
matrix. Let us for simplicity concentrate on the dependence of S on S' at a particular transverse 
position X. Suppose the wave function of the projectile at this transverse position has the form (we 
take there to be exactly one gluon at the point x) 


N^-l 


iTj,) = ^ Ca{x)\a,x) 


(2.18) 


a=l 


The scattering matrix operator when acting on T multiplies the gluon wave function by the matrix 

■S0- We thus have 

N^-l 

(T,|S|T,) = Y, Ca{x)C*,{x) S^’^ix) (2.19) 

a,6=1 

This generalizes easily to states with more than one gluon. Thus expanding the functional S in 
Taylor series in R we obtain the most general form compatible with its interpretation as the projectile 
averaged scattering matrix 


^[■«i = E E 

m=0 ni=l 


T-rm 




K., / 


( 2 . 20 ) 

where {a^} denotes the set a}, a^,..., , etc. Simply stated this is the 5-matrix of a state whose 

quantum mechanical amplitude to have rii gluons with color indices at the transverse 

coordinate Xi with i = I, ...,m is Cjai} {amj.(xi, ...,Xm)- Thus for example the coefficients of ’’diag¬ 
onal” terms with a* = 6* must be positive, since they have the meaning of probabilities to find the 
configuration with particular color indices in the projectile wave function. 
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Equations eq.( 2.15 , 2.17 , ^.20 ) determine the Hilbert space of the KLWMIJ Hamiltonian. We 
stress again that the norm in this Hilbert space as defined by eq.( |2.17 ) is very different from the 
standard quantum mechanical definition. The condition of positive probability density which in 
standard quantum mechanics is quadratic in the wave function, in our Hilbert space is instead a 
linear condition eq.( 2.20| ). This is a direct consequence of the fact that expectation values in our 
Hilbert space are defined as 


{0[p]) ^ / dpO[p] W[p\ 


( 2 . 21 ) 


which is linear in the ’’wave functional” W. 


The problem of calculating the S matrix as a function of rapidity in the present approach is 
posed in the following way. First find eigenfunctions of the KLWMIJ Hamiltonian, 


X 


KLWMIJ 


[R,5/5R]Gq[R] = UIqGq[R] 


The evolution of each eigenfunction is given by 

Gq{Y) = 


( 2 . 22 ) 


(2.23) 


Note that the eigenfunctions Gq[R\ do not have to satisfy eqs.( |2.17| , 2.20[ ). Instead we can take them 
to satisfy the standard Hilbert space normalization condition f DRG*q[R\Gq[R] = 1. Next expand 
at initial rapidity in this basis 




(2.24) 


The expansion coefficients are given in the standard form 


7 , = / DRGI[R]YP[R] 


(2.25) 


The expansion coefficients jq are of course such that is normalized according to eq.( |2.2[lD . 

The total cross section at rapidity Y is 


■s(r) = 7, 


(Y-Yo) 


I DSW^^[S{x)] G,[5(x)]. 


(2.26) 


Here we have assumed that the target weight functional W'^ depends only on S and have substituted 
the integration over a{x, x~) in eq.(^) by the integration over S{x) (ditto in eq. (|2.25D with respect to 
R and 5/5p). We note that the the unitarity of the scattering amplitude requires (barring miraculous 
cancellations) nonpositivity of all tOq. 


2.1 Correspondence with the standard RFT language. 

Before setting out to study the approximate eigenstates and eigenvalues of the KLWMIJ Hamiltonian 
we want to point out to direct parallels between our setup and the approach to Reggeon field theory 
utilized for example in |11]. 
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The standard RFT approach is based on a clear separation of s - and t channel with factorization 
into impact factors and t-channel exchanges. The t - channel exchanges are universal and do not 
depend on the nature of the projectile, while the impact factors are determined by the projectile 
wave function. This has direct parallels in our approach. The projectile averaged S - matrix eq.2. 2 ^) 
is determined entirely by the wave function of the incoming projectile. As explained above, each 
factor R in the expression eq.(2.20) corresponds to the scattering matrix of a single gluon in the 
projectile wave function. We will refer to these gluons as s - channel gluons for obvious reasons. On 
the other hand the eigenfunctions Gq[R] are completely independent of the projectile and are in this 
sense universal. They are determined solely by the evolution kernel ^ 

The eigenstates Gq[R\ are the exact analogs of, and in fact should be understood as, the t - 
channel exchanges. To see why this is the case, we note first that as is clear from eq. ( 2.201 ), Ti^[R] 
carries two sets of indices. The left index labels the quantum numbers of the incoming projectile 
state, and the right index labels the quantum numbers of the outgoing projectile state. These sets 
of indices include all conserved quantum numbers of the theory, in particular color and transverse 
momentum. In eq.(2.2C) we have indicated the color indices of each one of the incoming gluons a^, the 
total color representation of the incoming projectile state being determined by the product of these 
representations and analogously for the outgoing state. We assume for simplicity that the projectile 
state belongs to an irreducible representation of the color SUl{N) group and the outgoing state is 
in an irreducible representation of SUr{N). As for the ! transverse momentum, since the transverse 
coordinates of all the gluons do not change during the interaction, depends only on the difference 
of the initial and final transverse momenta k± — p±. On the other hand the Hamiltonian 
has all the same symmetries as discussed earlier in this section, and so the same quantum numbers 
are carried by its eigenstates Gq[R]. More precisely Gq[R] carry those quantum numbers which are 
not spontaneously broken by the vacuum of ^klwmij^ g^g below the SUl{N) ( 8 ) SU]i{N) 

group is broken spontaneously down to the vector subgroup SUv{N). Thus the set of quantum 
numbers q carried by Gq[R\ includes the representation of SUv{N), the transverse momentum and 
the discrete symmetry Z 2 0 Z 2 as discussed above. 

Referring to the calculation for the overlap coefficients 7 ^ in eq.( 2.25| ) we observe, that the 
integration measure DR is invariant under both SUv{N) and translations in the transverse space 
R{x) —> R{x + a). It is therefore obvious that the integral in eq.( 2.25| ) gives nonzero result only 
for those Gq which match the change in the appropriate quantum numbers of the projectile wave 
function in the scattering event. In particular g_L = k± — px and the SUv{N) representation of 
Gq is precisely that which is needed to change the SUl{N) representation of into the SU]i{N) 
representation. Thus the quantum numbers of Gq that contribute to eq.(2.24) are precisely those of 
the appropriate t - channel exchanges that contribute to the scattering. In this sense the eigenstates 
Gq[R\ indeed represent the t - channel exchanges with the target. 

Continuing this train of thought we also see that the expansion coefficients 7 ^ are analogous to 
impact factors of |^, which describe the coupling of the t - channel state with quantum numbers 
q to the s - channel projectile state. The analogy is indeed direct, however one has to distinguish 
between the impact factor of exact eigenstate of .^bich is Gq, and impact factors in the 

standard sense as they are used in []T|, namely the coupling of the projectile to a fixed number of 
t - channel gluons. To understand this distinction, we must realize that although an s - channel 
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gluon is identified with the factor R when it appears in of eq. (|2.20 ), the t - channel gluon is not 
identified with a factor R in Gq\R\ but rather with a factor 5/5p°' in its expansion. To see this we 
note that through eq.( 2.26| ) every factor of 6/6p°‘{x) in expansion of becomes a factor of the 

color field of the target a“(x) which is indeed naturally identified with the t - channel gluon. Thus 
the perturbative coupling of a fixed number n of t - channel gluons to a given projectile is given by 
the coefficient of the n-th order expansion of T,^[R] in powers of S/6p°‘{x). This coefficient, which is 
just the n gluon impact factor for the given projectile, can be written as 




(2.27) 


The latter equality follows from the explicit form of W[p] eq.( |2.15| ). 

This establishes the relation between the main elements of the standard RFT language and the 
basic quantities of our approach. To recapitulate: an s - channel gluon corresponds to a single factor 
R in at- channel exchange corresponds to an eigenstate Gq of the Hamiltonian 

a coupling of a given exchange to the projectile state is given by 7 ^; a single t - channel gluon is 
represented by a single factor bjbp'^ in the expansion of and an impact factor of a fixed number 
of t - channel gluons is given by Dq, eq.( 2.27|) . 

We also note that we can define the ’’rapidity dependent” function 


D^{Y)= Dpp^^{xi)...p^-ixn)WY[p] 


(2.28) 


which satishes the evolution equation 

d 


dY 


= / Dpp^^{xi)...p^ 


dp 


(2.29) 


This is the generalization of the 2 and 4-point functions considered in [T^ , |75[ |. The approach of 
jl^ , [T 9 I , usually considers this type of functions with the projectile consisting of a single dipole 
and customarily splits into the irreducible part and the reggeized part. To make contact with 
this approach we consider in Appendix A the evolution of for the projectile containing a single 
dipole in the leading order in 1/A^c- We demonstrate there that our approach in indeed equivalent 
to that discussed in |T^, ]T| and that the decomposition into irreducible and reggeized parts appears 


naturally in our approach as well. 


The rest of this paper is devoted to solution of the eigenvalue problem eq.( 2.22 ) in the approxima¬ 
tion where the KLWMIJ kernel is expanded around R= 1. This expansion can be formulated since 
i? = 1 is the classical minimum of the KLWMIJ kernel. Physically this expansion is best thought of 


as the approximation of small target. For a small target the overlap integral in eq.(2.26) is dominated 
by values of S close to unity. Thus the leading contribution at early stages of evolution comes from 
the wave functions Gq[R\ which are large at R close to unity. Those can be studied by expanding the 
KLWMIJ kernel around R = \ and keeping only the homogeneous term in the expansion in {R — 1). 
Note however that the condition that R is close to unity does not necessarily mean that we have to 
use the standard perturbation theory in ag, as this closeness does not have to be parametric in 
We discuss this point in Section 5. In the next section we will consider a simplified problem of the 
KLWMIJ kernel reduced to its dipole limit. 
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3. Reggeons in the dipole model 


A simplified version of the high energy evolution - the dipole model was introduced by Mueller in 
1^ . It describes the leading high energy behavior in the large Nc limit as long as the densities in 
the wave functions are not too large. As shown in |£7] the dipole evolution equation is obtained as a 
well dehned limit of the JIMWLK evolution. We can define the projectile dipole ’’creation operator” 
s and the target dipole ’’creation operator” r: 


s{x,y) = ■^Tr[Sl,{x)SF{y)]-, 


1 


r{x,y) = -^TrlR^p^ix) Rpiy)] 


(3.1) 


where F indicates the fundamental representation. The term ’’creation operator” is not mathemati¬ 
cally perfect - it would be more appropriate to call it the ’’dipole field”, since it is Hermitian. We will 
nevertheless keep to the tradition of calling it the creation operator in this paper. The ’’annihilation 
operators” in this parlance are ^ and ^ for projectile and target respectively. 

If the target weight function is a function of s only, that is = IT^[s], the action of the 
JIMWLK kernel on it in the large Nc limit is equivalent to the action of the dipole kernel (see 
also 1 ^ ): 


X 


JIMWLK 


W^[s\ = 




W^\s\ 


(3.2) 


The dipole kernel was found in |64, 65] by reformulating the original Muller's model and was later 
obtained directly from eq.(|3.2|) in j^j 


dipole 

Xs 


S, 


Ss 


Mx,y,z [-s{x,y) + s(x,2;)s(y,2;)] 


'x,y,z 


5s{x,y) 


with M - the usual dipole emission probability 


, Ols 
Mx,y,z = ^ 


(x - yf 


2 7r (x — z)2 {y — zY 


(3.3) 


(3.4) 


Eq.(3.2) is derived under the same assumption as the JIMWLK equation - namely that the target 
is dense. The dual form of this kernel describes the evolution of a hadronic state (be it a projectile 
or a target) which contains a small number of dipoles in its wave function. This form is derived from 
the KLWMIJ kernel assuming that r is the only degree of freedom in the wave function (see also 


X 


KLWMIJ 


IT[r] = 


’ 5r 


W[r] 


(3.5) 


The calculation of the cross section of a projectile on a target which both are made entirely of dipoles 
is given in analogy with eq.(2.26) by 


5z,(y) = 


(Y-Yo) 


Ds IKyJs(x,y)] gq[s{x,y)]. 


(3.6) 


where the wave functions Qq satisfy 


=uqgq[r] 


(3.7) 
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We start by expanding the dipole Hamiltonian around the classical solution r = 1. Let us denote 
the dipole creation an annihilation operators by^ 


= 1 - 


r ; 


^ 6r 


The dipole Hamiltonian is 


^dipole ^ ^ 

The free Hamiltonian Hq is quadratic in the creation and annihilation operators 


Ha = 


Jx,y,z 


-S{x,y) + d\x,z) + d){z,y) d{x,y) 


The Hamiltonian of the interaction Hj is given by 


Hi = - M^^y^zd){x,z) d) {y,z) d{x,y) 

Jx,y,z 


(3.8) 


(3.9) 


(3.10) 


We now consider the spectrum of the free Hamiltonian Hq. There are of course no new results in 
this calculation. It only serves to establish in a very simple setting the correspondence between the 
Reggeon field theory and the KLWMIJ/JIMWLK approach that we want to make explicit in this 
paper. This is an easy exercise, since the Hamiltonian is very similar to harmonic oscillator. The 
vacuum of Hq is |0) such that d|0) = 0. The wave functional of this state is simply a constant which 
does not depend on r. The first exited state contains one dipole and is a superposition of the basis 
states d'l^(x,y)|0). We take it in the form 


9q = 


J 

J Xd 


1 


■'ipq{x,y)d^{x,y) |0) 


>x,y {x - y)^ '' 

Acting on gq by Hq we find that the wave function should satisfy the BFKL equation 
/ M^,y,z [-^|Jqix,y) + 'tpq{x,z) + 'tpqiy, z)] = iOqljjq{x,y) 


(3.11) 


(3.12) 


The eigenfunctions should vanish when two transverse coordinates coincide, since the dipole of zero 
size is indistinguishable from vacuum. Solutions of the BFKL equation which satisfy this condition, 
the so called Moebius invariant solutions, are well known. They are the eigenfunctions of the Casimir 
operators of conformal algebra |p^ 


V’g(x,y) = E^'''{x - p,y - p) = 


(x - y) 


{x - p){y- p) 


1 + n 


+iu 


{x - y)* 


{x - pY {y - pY 


-+iu 


(3.13) 


where the complex coordinate x is dehned as x = xi + 1 x 2 The index q denotes the conformal spin 
n, v, as well as the degeneracy vector p. The eigenvalues are 


xiq = 2asx{n,i^), 


(3.14) 


^Although we are abusing notation slightly by calling both r and 1 — r the dipole creation operator, we hope this 
does not cause confusion. 
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where 


(3.15) 


X{n,u) = V’(l) -Irtp 


1 + |n| 


+ iv 


--iP 
2 ^ 


1 + |n| 


— lU 


These wave functions satisfy the completeness relation |68| 


°° roo !■ 16 

{2-nf 5{x - x) 5{y - y) ='^ dv d^p-—^^-^-——4^ E"”^{x - p,y - p) E^’''*{x - p,y - p) 

J-oc J \x y\ \x y\ 

(3.16) 

By diagonalizing Hq in the one dipole Hilbert space we have found such a linear combination of s - 
channel dipole states which couples precisely to a given ’’state” (fixed (n, u)) on the BFKL Pomeron 
trajectory. Note that this does not mean that the coupling of this state to the projectile and the 
target is via a two gluon exchange only. Since the wave functional pq depends linearly on r, and r has 
all order expansion in powers of S/6p, the coupling of pq to the target is the full eikonal coupling which 
includes multiple scatterings. Nevertheless this t - channel exchange with all its multiple scatterings 
evolves exponentially in rapidity with a fixed exponent 2dsx(^) ^)- 


The exited states at the next level contain two ’’reggeized dipoles”. Since in the dipole model all 
dipoles are treated as independent degrees of freedom and [d^{x,y),d{u,v)] = —6‘^{x — u)6‘^{y — v), 
we immediately infer the two dipole spectrum 


ddoPqi,q2 ~ [^qi ‘^q2\9q\,q2 


(3.17) 


where 


9qi,q2 — 


' X,y,UyV 


(x — yY {u — vY 


Yqi {x, y) Yq 2 {u, v) d^ {x, y) d^ (n, u) |0) 


(3.18) 


These eigenstates correspond to double Pomeron exchange. This corresponds to the leading 
order in BKP state (which couples to a projectile consisting of two dipoles) which contains four 
reggeized gluons in the t-channel. Note however that as discussed in section 2, this exchange does not 
contain a fixed number of t-channel gluons, but rather includes all multiple scatterings. If we chose 
the projectile S to be equal to Pqi,q 2 , such a projectile state in the present approximation would not 
interact via a single Pomeron exchange, but only via the double Pomeron exchange. As discussed 
at the end of this section however, such a projectile is not physical, as its wave function contains 
negative probabilities. 

The generalization to the multi-dipole states is straightforward. The eigenstates and eigenvalues 
of Hq are 

HoPqi,...,qn — {^qi E ■■■ E tOqY Pqx,...,q„ (3.19) 


with 


9qi,...,qn 


f Yqiixi,yi) 

'{xiUyi} - yi)^ ' 


Yq-n {Xn-i Vn) 

{Xn UnY 


dYxi,yi)...dYxn,yn)\0) 


(3.20) 


Thus the dipole limit of the KLWMIJ hamiltonian reproduces the leading large Nc spectrum of 
the noninteracting multipomeron exchanges. Note that the eigenstates contain all the nonforward 
BFKL amplitudes and thus scattering amplitudes at arbitrary transverse momentum transfer. 
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This of course is simply a recasting of well known results in our framework. It is therefore 
also obvious that the expansion around r = 1 violates unitarity. While the complete dipole model 
Hamiltonian is unitary and when solved nonperturbatively must lead to unitary S - matrix^, the 
Hamiltonian Hq does not. One manifestation of this is the fact that the BFKL trajectory has an 
intercept greater than one, and therefore the scattering amplitudes grow without bound. The same 
fact also manifests itself in a somewhat different fashion. In particular the eigenfunctions that we 
found are all (except for the vacuum) zero norm states with respect to the norm defined in eq.( |2.17| ). 
All the wave functions are homogeneous polynomials of 1 — r, and thus they all, except for the 
constant vacuum wave function, vanish at r = 1. 


Put in other words, when interpreted as dipole states in the s-channel, Qq contain negative 
probabilities. The condition eq.( 2.2C| ) can be adapted to the dipole limit in a straightforward manner. 
It states that in the expansion of S[r] in powers of r(x, y), individual coefficients represent probability 
densities of finding a particular dipole configuration, and thus all coefficients of r(xi, yi)...r(x„, y^) 
must be positive. On the other hand consider the first excited state wave function 


99 = [ [rix,y) - 1] . (3.21) 

j x^y y) 

The coefficient function ijjq as defined in eq. (|3.13|) is not real. Due to the degeneracy of the eigenvalues 
eq.(ra), n —> —n and v —v one can construct real combinations of these functions 
and i{'ipn,u — '4’-n,-u)- These combinations however are not strictly positive. Thus if we were to take 
S[r] = gq[r], this wave function would have negative probabilities to find the dipole at some points 
(x, y). The only eigenfunction which is positive everywhere is the one that corresponds to the forward 
scattering BFKL amplitude, n = v = 0. Even this one however contains negative probabilities. Since 
go is proportional to r — 1 and not r, it has positive probability to contain one dipole (at any position 
(x,y)) but an overall negative probability to contain no dipoles at all. Similarly, the second excited 
state has a positive probability to contain two dipoles and no dipoles, but a negative probability to 
contain a single dipole. 


This curious behavior is not by itself a problem, since one can construct a positive norm state 
with all positive probabilities at the initial rapidity by taking appropriate superposition of these zero 
norm states and the vacuum state. This is also the reason why we have to keep all the eigenfunctions 
that we have found in the preceding analysis and not just the one with re = = 0 - we need to 

superpose all of them in order to construct an arbitrary positive definite scattering matrix S[r]. The 
problem arises however at later times (higher rapidities) since the BFKL intercept is greater than one. 
The S-matrix at later times will be dominated entirely by the component in the superposition which 
has the largest eigenvalue, and this state by itself contains negative probabilities. The linearized 
evolution therefore preserves the overall norm of the state (since S[l] does not evolve with rapidity), 
but generates probabilities which are greater than unity as well as negative probabilities. Thus unless 
the unitarizing corrections due to the triple Pomeron vertex are taken into account, the evolution 
violates unitarity not just by driving the value of the total scattering probability above unity, but 
also by making the probabilistic interpretation in the s - channel impossible. 

'^Within additional approximation of no target correlations (see for alternatives), finding a non-pertnrbative 
solution of the dipole Hamiltonian is equivalent to solving the non-linear BK equation. This question was addressed in 
many papers ©US §,g. The result is consistent with the unitarity constraint. 
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We certainly believe that this is an artifact of the perturbative approach and that the interaction 
will cure this problem when consistently taken into account. It is not our aim in this paper to solve 
nonperturbatively the KLWMIJ hamiltonian (we wish we could!) but rather to relate the perturbative 
approach to the standard techniques which also operate with nonunitary amplitudes. We therefore 
will not have much to say about the effects of the interaction. We only note that the first perturbative 
correction to the single Pomeron state can be obtained by acting on the state by the perturbation 
Hamiltonian Hj: 

Higq = [ ^pq{x,y)d\x,z) d^{y,z) |0). (3.22) 

Jx,y,z I® y) 

This state can be projected onto a two-particle state gqi,q 2 - The resulting matrix element is 


{9c 


gi,<i2 


Hi \gq) = - [ 

Jx,y,z 


M, 


x,y,z 


(x - y)^ (x - z)-^ {y - zY 


Yq{x,y)Yl (x,z) 


(3.23) 


Eq. ( 3.23| ) is the well known result for the first unitarity correction to the BFKL Pomeron due to 
the QCD triple Pomeron vertex. In the dipole model and beyond this has been studied in detail in 
Refs. |I1, 12, [^, 1^, 28, 23]. It is a satisfactory feature of the present approach that this result 

is reproduced automatically. 


We therefore conclude that the dipole model Hamiltonian considered as a second quantized 
2-|-l dimensional field theory describes in the perturbative expansion the large W BFKL Pomeron 
with the self-interaction given by the triple pomeron vertex. It is tempting to think about a given 
eigenstate of Hq as a ’’projectile” with specific conhguration of dipoles in the s-channel chosen so 
that it interacts with any target only by the specific n-Pomeron exchange. The wave function of such 
a projectile however contains negative probabilities. It is thus more correct to think about a given 
eigenstate as of a t - channel exchange, as discussed in Section 2. Any physical ’’projectile” with 
all positive probabilities can be constructed at initial rapidity as a superposition of these eigenstates. 
The evolution generated by Hq however generates negative probabilities at large enough rapidity. 


In the next section we perform a similar analysis of the full KLWMIJ Hamiltonian beyond the 
dipole limit. 


4. Reggeons in QCD 


4.1 The Hamiltonian 


We start with the KLWMIJ Hamiltonian written in the form eq.( p.l4 ). It has a classical minimum 
at R = 1 and in the rest of this section we consider expansion around this configuration. Introducing 
R = 1 — R we write 

(4.1) 


KLWMIJ _ TT , TT 

X — Ro + Ri 


The “free” Hamiltonian Ho is the part homogeneous in R —, and it naturally splits into the real and 


virtual parts, the latter appearing due to the normal ordering in eq.(2.13) 


Ho = H§ + hX 


(4.2) 
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The real part is 


— Kx^y^z 


f 

d n , ~ ~ ^ 


~ d 

\ 

6 n , ~ ~ N 

< — 2tr 

-^T^{Rx - Rz) 

16 rI \ 

tr 

{Ry - Rz)T^^ 

6R^ 

> + tr 

-^T^Rx - Rz) 


0 tr 


6R 


r“ {Ry - R, 


+ tr 


{Rx - Rz)T^ 


5rI, 


tr 


(Ry - Rz)T^ 


5Rl 


The virtual part is 


< = -2 


Kx 


— tr 


T’^iRx - Rz)T^ 


6Rl 


+ N tr 


{Rx - Rz 


The interaction Hamiltonian has the form 

5 


Hi = -K 


x^y,z 


tr 


6R. 


T“ {Rx - Rz) 


rI^ tr 


{Ry - Rz)T^ 


6Rl 


5rI 


(4.3) 


(4.4) 


+ 


I 


Kx,x,z [Rf + &T 


tr 


T^{Rx - Rz)T^ 


5RI. 


(4.5) 


Here we have used the unitarity of the matrix R to rewrite the first term in the form which looks 
like a 2 —> 4 vertex. We have dropped the normal ordering sign in the above formulae for notational 
simplicity, but all the functional derivatives are understood not to act on the fields R in the kernels 
themselves. 


We view the Hamiltonian eqs.(4.1.4.3.4.4,0|) as defining the fully second quantized Reggeon 
Field Theory. In the following we will refer to it as either KLWMIJ or RFT Hamiltonian. 

As noted in Section 2, in principle the separation of the full RFT Hamiltonian into a free term 
and an interaction is not unique due to the unitarity of the matrix R. Our guide is to define this 
split in such a way that Hq and Hj are separately ultraviolet and infrared safe. We discuss this point 
further in Section 6. 

Analogously to Section 3, we can think of R and 5/5W as s - channel gluon creation and 
annihilation operators respectively. Consistently with this terminology we will call the approximation 
based on the Hamiltonian Hq in eq.(^) - partonic approximation to RFT. In this approximation 
the s - channel gluons emitted throughout the evolution scatter independently and the scattering 
amplitude of a gluon configuration is taken to be the sum of scattering amplitudes of the individual 
gluons. This is encoded in the Hamiltonian Hq since it is homogeneous in the gluon creation operator 
and thus preserves the number of s - channel gluons that contribute to the scattering amplitude 
throughout the evolution. We will discuss the physics of this approximation as well as its relation to 
the standard BFKL/BKP approach in the next section. 

The vacuum of the free Hamiltonian Hq is clearly the state |0) such that 




| 0 ) 


= 0 


(4.6) 


This is the state with no gluons and corresponds to a trivial ’’empty” projectile. The higher states 
are obtained from the vacuum by successive action of the creation operator R. Although Hq is 
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homogeneous in the creation and annihilation operators, it is significantly more complicated than 
the corresponding Hamiltonian in the dipole model. It contains not only terms of the type a^a, but 
also a^a^aa, and its spectrum is therefore not simply that of a harmonic oscillator. 

In the diagonalization of Hq we will hnd very helpful to use its symmetries. The expansion 
around R = 1 breaks the SUl{N) 0 SUii{N) symmetry of the full KLWMIJ kernel down to the 
diagonal SUv{N): 

R{x) R{x) U . (4.7) 

Thus the eigenstates of Hq can be classihed according to the representations of this SUv{N) group. 
As discussed above, and as we will see explicitly below, the representation of SUv{N) corresponds 
to the color quantum number of the t - channel exchange. 

Another interesting symmetry of Hq and Hi is the discrete symmetry 

R{x)^R^x) (4.8) 

This transformation interchanges the row and column indices of the matrix R. Physically the first 
index of the matrix R corresponds to the color index of the incoming s - channel gluon, while the 
second index to the color of the outgoing s - channel gluon. Thus the interchange of these indices 
interchanges the s and u channels. The discrete symmetry eq.( |4.8D can be therefore identified with 
the signature discussed in the context of the Reggeon theory®. The eigenstates of Hq fall into even 
and odd states under this transformation. 

Finally another discrete symmetry of the KLWMIJ Hamiltonian is charge conjugation (C-parity). 
In QCD the action of the charge conjugation symmetry is conventionally defined as 






(4.9) 


where r“ are the generators of the fundamental representation. This translates into the following 
transformation of the fundamental eikonal factor 


Rf{x) R*f{x) 

while the adjoint eikonal factor transforms as 


(4.10) 


(4.11) 


Here the matrix C is the diagonal matrix which can be conveniently written as 

= 2tr(r“r*®) 


(4.12) 


The transformation eq.(4.11) flips the sign of the matrix elements of R in rows and columns whose 


indices correspond to the imaginary fundamental generators r“. This symmetry is also unbroken by 
the configuration R = 1. 


®We dub “signature” the behavior of the amplitude under the s — u crossing. This is not quite the standard definition 
of the signature factor as used in the partial wave analysis, but coincides with it for scalar exchanges. We thank Genya 
Levin for pointing this to us. 
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4.2 One-particle state: reggeized gluon and other Reggeons 


In this subsection we discuss the ’’one particle states” - the eigenstates of Hq linear in the s- channel 
gluon creation operator R. We will use the terms ’’one particle state”, ’’two particle state” etc. in 
this sense throughout the rest of the paper. 

The one particle eigenstates represent t - channel exchanges which in the partonic approximation 
couple directly to the one s - channel gluon states. We are looking for one-particle eigenstates of the 
form 

(4.13) 

J X 

The states in this sector are multiplets of the SUv{N) symmetry, and we use cumulative index A to 
label both the representation and the particular state vector inside the representation. The tensors 
essentially are projectors from the product of two adjoints to the representation A, see below. 

The eigenvalue problem we are to solve is 


= /3,(A)G^ 


(4.14) 


First we compute the action of Hq on R 

jK^,x,z [{T-T^)cd{Rf - Rf) - N{Rf - 


Tj fycd _ zjV _ o 

iif) ri^ — iiQ — Z 


)} 


(4.15) 


The color structure can be easily diagonalized. To do that we need solutions of (no summation over 
A is implied) 

(4.16) 

This is solved using the decomposition of the product of two generators into the projectors of SU(N) 
representation (see eq.( C.22 ) of the Appendix C) 

(r“r'),, = N = N (^P[l] + -P[8a] + -P[8s] - -P[27] + P[R 7 ]j ^ , (4.17) 

The explicit form of the projectors P is given in the Appendix C. We follow in this paper the notations 
and nomenclature of [^]. Thus the different representations are labeled by the dimensionality of 
their counterparts in SU{3), except for TZj which does not exist in the SU{3) case. Clearly choosing 
r] = P[A] solves equation ( |4.16 ) for any representation A that can be constructed from the product 
of two adjoints. Thus we have seven eigenfunctions corresponding to 1, 8 a, 8s, 10, 10, 27, and TZy 
representations of SUv{N) with the eigenvalues 


Ai — 1, 


■^8a - 


1 


^Ss - 


1 


Alo — — 0 , A27 — — 


1 

iV’ 


1 

N' 


XTlr = iT- ( 4 . 18 ) 


Note that each projector P[X\ corresponds to a t-channel exchange in the color representation 
A. This interpretation follows directly from the fact that R is the T-matrix of the s-channel gluon. 
For ^[A] the color indices of the incoming and outgoing gluon are projected onto the representation 
A. Physically this can happen only due to an exchange with the target by an (t - channel exchange) 
object in this representation (Fig. |l|). 
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Figure 1: t-channel projection of Wilson lines. 


Once the color structure is diagonalized, the equation for follows: 


- 2N{1 - A) / = /3,(A)^,(y) 

J Z 

Substituting 'h in the form of a plane wave 'l'q(y) = we get 


(4.19) 


a.o N 


TT^ 


(1 - A) / (fk 


I 

J u 


{q — k)"^ 




(4.20) 


where /r is the infrared cutoff. The momentum q is the momentum transfer in the scattering process. 
For the eigenvalue (3q we obtain (cis = ag N/n) 


(Xc 


/3g(A) = - —(1 - A) d^k- 

TT ' 


- Os (1 - A) In ■ 


(4.21) 


k‘^{q-kY ■ 

In the last equality we have assumed q ^ fj,. For q ^ 0 the eigenvalue vanishes quadratically with q: 

..2 


Q 

/3q(A) >0 CTs (1 A) 2 


IJ- 


(4.22) 


Note that the coefficient 1 — A is simply proportional to the second Casimir of the given representation 
(tabl. 1^: 

2iV(l-A) = C2. (4.23) 

Eq. ( [4.21| ) provides an expression for the trajectory of the Reggeon in the channel with given color 
exchange |11[| (see also ||^). The TZ = 8a channel is special since it has the quantum number of the 
gluon. Indeed we have 

PqiSA] = - 


asN q^ 

- In —^ 

2 vr iJ? 


which is the standard reggeized gluon trajectory. We will refer to it as the reggeized gluon or the 
/-Reggeon. Its signature is negative, since the projector P[8aT^ is antisymmetric under the exchange 
of a and b. The /-Reggeon has a positive signature ’’brother” in the TZ = 8s channel, which we will 
refer to as the d-Reggeon. Since Asg = Ag^, its trajectory is degenerate with that of the reggeized 
gluon, / 3 g[ 85 ] = Pq[8A]- The corresponding wavefunctions are 


C^^Aab ^ f ^iqx jabk jkcd 

J X 


f lP<^{x) 

J X 
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Thus in our approach the reggeized gluon is naturally identified with an eigenstate whose eigen¬ 
function is linear in the Wilson line R projected onto the antisymmetric octet representation in the 
t-channel. Note however that the identification is not operatorial - we do not have an operator in 
the Hilbert space which we could identify with Reggeized gluon. The correspondence is rather on 
the level of the (perturbative) eigenstate of the KLWMIJ Hamiltonian. 

Note that the singlet state in the t-channel {TZ = 1) does not contribute to the evolution of the 
cross section. The scattering amplitude due to the exchange of this ”s-Reggeon” is flat with energy 
since Ai = 1. Note also that the Reggeons corresponding to 10 and 10 exchanges are degenerate. 
Thus if the coefficient of the projector in the expansion of 5][R] is equal to that of at initial 
rapidity, it remains so throughout the evolution. Nevertheless the 10 and 10 Reggeons appear as 
distinct eigenstates of the RFT Hamiltonian. 


The signature of the 10 and 10 Reggeons as well as of the reggeized gluon is negative since the 
respective projectors are odd under the exchange of one pair of indices. The other three Reggeons - 
85 , 27 and 1Z^ have positive signature. Different types of Reggeons (including the singlet) and their 
relation to signature were discussed in 0, MM- 

The eigenvalue at —> 0 for all the Reggeons vanishes. For the reggeized gluon in the present 
approach this is the consequence of the Goldstone theorem. Our expansion breaks spontaneously 
the SUl{N) (g) SUji{N) symmetry down to SUv{N). As in any other quantum field theory, in the 
RFT such breaking requires the appearance of — 1 signature odd Goldstone bosons in the adjoint 
representation of SUv{N). These are precisely the quantum numbers of the reggeized gluons. The 
fact that all the other Reggeons also have vanishing eigenvalue at = 0 from this point of view is 
surprising. It may signal the presence of a larger broken symmetry group which we were not able to 
identify. 


The effect of the interaction Hamiltonian Hj is to mix the one and two Reggeon states. In 
particular acting by the perturbation Hj on the one Reggeon state sector we obtain: 


HiG^g = [ 

J X.: 


K^,x,zRf T^{R,-R,)T^ 


- ab 


^q{x) 


(4.24) 


Projecting the above state onto the analog of the noninteracting two-particle state using the projector 


' H,Q2 — 


L 




6R^^{x) 






_ 

5R'^\x) 


(y) 


we obtain 

= vll qi, q2) (4.25) 

with the Reggeon ’’splitting vertex” given by 


Vig^2g{q]qi,q2) 


I 

Jx,y 

2^ 


Ay A{x)%^{y) [4^,(x) + 4/,(y)] 
5^(9 - qi - q2) 


In4+ln4' 


(4.26) 


In accordance with the signature conservation not all of the transition are allowed. For example the 
/-Reggeon mixes only with the / and d states. 
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4.3 Two-particle state 


Our next task is to find eigenstate in the two particle sector. We will limit ourselves to the SUv{N) 
singlet states which correspond to color singlet exchanges in the t-channel. It is convenient to label 
these states Gq^ by their eigenvalues with respect to the discrete symmetries: signature S, charge 
conjugation C and space parity P which for two particle states is equivalent to the symmetry of the 
eigenfunction with respect to transverse coordinate exchange. The index q cumulatively denotes the 
rest of the quantum numbers. We will refer to the C-parity even eigenstates as Pomerons, while to 
the C-parity odd eigenstates as Odderons. The most general Pomeron state has the form 

6 

= R^Hu)P‘^{v)^liu,v)\oy, (4.27) 

i=i 

Here the projectors P* = P[JZi] are the same as in the previous subsection and the index i runs over 
representations 1, 8 a, 8 ^, 10 -|- 10, 27, TZj. We have defined the projector onto the representation 
10 -|- 10 as -|- R^^]. Note that the combination — P^^ does not appear in the Pomeron sector, 
since it has negative C-parity. 

In the Odderon sector we have three possible sets of eigenstates which all differ in at least one 
quantum number: 

C" = Z+“ f R^\u) R^^{v) 4>- (n, u) |0) ; 

R U,V 

C+ = Z-“ [ ip\u) P\v) ch+(n, v) |0) ; 

J u,v 

Gy = [P^^ - P'^ti f |0) ; (4.28) 

J u,v 


a 


b 


s 




c 


d 


a 


b 



a 


b 



c 


d 


a 


b 


c 


d 


Figure 2: s-channel projection of two Wilson lines. 


The additional tensors are defined as (see Appendix C). 

= 2 (^fbakdkcd ± dbakfkcd^ (4.29) 

The projectors P* are symmetric under the exchange of indices P^jf^ = P^cjf while the tensors 
are antisymmetric under the exchange of both pairs of indices. Therefore the ^l{u,v) and 4'“(u, u) 
are functions symmetric under the interchange of u and v {P = +1) while ‘h^(u, v) are antisymmetric 
functions (P = — !)• 
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Eq. (4.27) and (4.28) have a simple interpretation from the s-channel point of view (Fig. 


^ ^ ^ ^ ^ _ ___ In 

two incoming gluons with indices ac are projected onto representation TV'. The two outgoing 
gluons (with indices bd) are projected onto the same representation TZ^. Thus the color representation 
of the two gluon state remains unchanged during the propagation through the target. Additionally 
the invariant tensor P* projects the product of the incoming and the outgoing representations TV 
onto a color singlet. Thus these amplitudes correspond to a color singlet exchange in the t-channel. 
It can be easily checked that these eigenstates are even (C = +1) under the charge conjugation 


transformation eq.( 4.11 ). The tensors P* are symmetric under the interchange of the incoming and 
outgoing indices; P*ac = P^^d- Hence all are positive signature eigenstates {S = +1). 

In \ the incoming and outgoing representations are both adjoints, but they have different 

charge conjugation properties. These eigenstates correspond to the change of 8 a into 85 (or vice versa) 
during the interaction with the target. Again one can straightforwardly check that both G^( ) 

eigenstates are odd under the charge conjugation eq.( 4.11 ). They are therefore both associated 
with the Odderon exchange. The properties of these two sets of eigenstates under the signature 
transformation are however different. The tensor is odd under the interchange of the incoming 
and outgoing indices and therefore the set G~^ ^ contains negative signature eigenstates®. These 

correspond to the standard signature odd odderon exchange. On the other hand the tensor Z~ 
is even under the interchange of the incoming and outgoing indices. The eigenstates G^ ^ ^ thus 

correspond to a signature even Odderon exchange. 

The set of eigenstates G~ ^ is a peculiar one. These states have negative signature and 
negative charge conjugation but are described by a symmetric wavefunction. Because the Hamiltonian 
preserves parity, charge conjugation and signature, the four sets of eigenstates of eqs. ( 4.27 , 4.28 ) do 
not mix. 

The s-channel representation eqs. ( |4. 27 , 4.28 ) is natural for discussing the scattering process from 
the s channel point of view. It emphasis the color structure of the gluon ’’dipole” state which 
propagates in the s - channel. It is however also instructive to consider an alternative t-channel 
representation (Fig. ^). Instead of projecting the color indices of the two incoming gluons into a 
definite color representation we can project the incoming and outgoing indices of the same gluon. 
This corresponds to writing the wave function of the eigenstates in the form 


Q+ (++) _ \ ' pi cd 


TP\u) R‘^\v)^\{u,v) 


ycd 


i 


ab 


(4.30) 


In this expression a given term in the sum eq.( 4.30 ) represents the process whereby each one of 
the s-channel gluon exchanges color representation TV with the target, and additionally the two 
77.* representations in the t-channel form an overall color singlet. The two expression eq(4.27) and 


eq. (Pa) are simply related by the change of basis, since either set of projectors (s-channel or t- 
channel) form a complete set. The t-channel coefficient functions 'kj are linearly related to the 


functions 'k*: 




= 




\j/^ 

I ^ t ^ 


(4.31) 


®We hope that our convention in which the tensor defines the signature odd state is not hopelessly confusing. 
Although this convention is somewhat unfortunate, we have decided to stick literally to the notations of |^] rather 
than risk even greater confusion. 
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Figure 3: The t-channel projection of two Wilson lines. Each one of the two s - channel gluons exchanges a 
given color representation with the target. Additionally the two exchanges form an overall color singlet. 


where the crossing matrix C is given in the Appendix C. 

Analogons transformation can be performed for the Gq ^ \ \ Gq ^ states. The s 

channel Z~^ tensor remains the tensor in the t-channel, while the s-channel Z~ tensor becomes 
in the t-channel the difference of the projectors — P^°]. The s-channel [P^*^ — P^^] becomes Z~ 
in the t - channel. 

S iCP) 

This t-channel representation of the fnnction Gq^ Ms snggestive of the interpretation of the 
exchanged states in the t-channel as bound states of Reggeons discussed in the previous subsection. 
In fact this turns out to be a very convenient interpretation since as we will see the eigenstates of Hq 
correspond to a fixed t-channel projector P* (that is only one y) is nonvanishing for a particular 
eigenstate). 

The eigenvalue problem in the two particle sector is 


(4.32) 

Our goal now is to reformulate Eq. ( 4.32|) as equations for and <h. This is achieved by applying 
to Eq. ( 4.32 ) of the operators 

<5 <5 


pi _ 

(5PM(x) (5PP(?/) 


and Z'"^^ - ^ -- 

^ 5PM(a;) (5PP(y) 


This leads to the following equations 

6 6 


pj _ 

(jpM(a;) 5PP(y) 


PoG+(++) = 2u;qtr[P^]n{x,y) = 2u;q A>i{x,y) (4.33) 


and 


--4- HoGj^—'^ = ^4a;„(Ai2 _ 4:)Ds<^^{x,y) 

“^<5PM(x) 5PP(y) ) s K 

Z~ ^4--4- Ho G- H+) = 4a;„ {N^ - 4) Pg K (x, v) 

“^5PM(x) 5PP(y) ^ t a sK ,yj 


(4.34) 


(4.35) 


where Dj = dim\R,j\ is the dimension of the corresponding representation. The eigenvalues Uq 
and the functions <I> are found by solving the system of homogeneous equations Eq. ( 4.33 ) and 
Eq. ([4.34}) . The details of the calculation of the left hand sides of Eq. ( |4.33D , Eq. ( |4.34D and Eq. ( 4.35|) 
are given in the Appendix B. The results are summarized and discussed in the next two subsections. 
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4.3.1 Pomeron and family. 


Let us consider first the equations for the C-parity even sector. Equations eq.( 4.35 ) can be written 
as (see Appendix B) 




2 Arx,y ,2 - Kz^y,x'^l{y,z) +2 6{x-y) / Kz,u,x%{z,u) 


U) 


+ E / [‘^Kx,x,.-^l{x,y) - K,,,^y¥,{x,z) - K,,,,x'fliy,z)] = ^^iix,y) 

(4.36) 

where Sj = 0 if f is a symmetric representation, Sj = 1 if i is antisymmetric, and the matrix C is 
given in eq.( p.25 ). No summation over j is implied. Using the relations eq.(4.31) one can show that 
this system is diagonal in the t-channel 


L 


‘^Kx,y,z '^\{x,y) - Kx;,z,y^\{x.,z) - Kz,y,x^liy,z) + 2 6{x-y) / Kz,u,x^l{z.,u) 


- / [2Kx,x,z'^l{x,y) - Kz^z,y'^l{x,z) - Kz:,z:,x'ifl{y,z)] = 


iO 


mx,y) (4.37) 


2iV(l - AO 

with Aj given by eg.( 4.18 ) and A^Q_|_j^ = 0. Eq. ([4.37 ) is the non-forward BEKL equation in the 
coordinate representation. In the Eourier space 


0(^1, ^ 2 ) = / mx,y)e 

Jx,y 


ikix-\-ik2y 


the equations eq. ([4.37|) read 


U) 


[2(1 - Ai) 


- (5{ki) - Pik2) 


0(^1, ^ 2 ) = / K2^2iki,k2-,qi,q2)'il^{qi,q2) (4.38) 


' 91,92 


The kernel K 2 ^ 2 {ki, k 2 ',qi,q 2 ) is real part of the of the standard momentum space BEKL kernel 


a 


K2^2{ki,k2]quq2) = - — 


{qi - ki)i gi, 


{q2 — ^ 2)1 _ ^ 
2tx [{qi - kiY ql\ [(^2 - ^ 2 )^ ql 


6 ‘^{ki + k2 - qi - q2) 


(4.39) 

The BEKL trajectory is given by Eq. ( p.l4 ). Thus we find five independent singlet exchanges 

in the t-channel with the eigenvalues 


J = 2{l - (4.40) 

They have a natural interpretation as ’’bound states” of the Reggeons discussed earlier. The com¬ 
pound state of two reggeized gluons is the standard BEKL Pomeron with 

oo[Sa.Sa] = (4.41) 
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In addition we have the following nontrivial bound states (Fig. Q) 


w[ 85 , 8 s] = , a;[10 + 10, lO + TO] = , 

a;[7^7,7^7] = 2 (^1 - a;[27,27] = 2 (^1 + . (4.42) 


y-:- 

1 Pomeron 



8 a 8a 
w = ff> bfkl 





Pomeron 





8s 8s 

0) = CO bfkl 













27 27 



(0 — 2 ( 1 + 1/N ) CO BPJ^L 


Figure 4: t-channel bound states of reggeized gluons. 

The Pomeron that would correspond to exchange of a pair of singlet representations has zero 
eigenvalue just like its ’’constituent” s - Reggeon. Thus the contribution due to this exchange to the 
cross section is flat with energy. 

Note that the eigenvalues in the 10 + 10, 27, and TZ-/ channels are proportional to the second 
Casimir of the respective representations and are therefore greater than that of the BFKL Pomeron. 
The reason is that these exchanges require at least four gluon coupling in the t-channel as opposed 
to the BFKL Pomeron which in the Born approximation consists of two t-channel gluons. Thus the 
coupling of these Pomerons in perturbation theory is of order and they have the same nature as 
the two BFKL Pomeron exchange. We will discuss this in more detail in the next subsection. 

It is interesting to note that the eigenvalue of the bound state of two 27-Reggeons is higher 
than twice the BFKL value. Although four t channel gluon states with intercept greater than two 
Pomerons have been discussed in the literature [0, 0], the correction we find is higher. It is of order 
l/A’c and not 1/N^ as may be naively expected. The reason for this is that in the ^ oo limit the 
27 and 77-7 Reggeons are degenerate. In this situation the finite correction naturally starts with 
order 1/Nc as the 1/Nc perturbation theory is degenerate. This is also the reason why the 0(1/A^c) 
corrections to the TZy and 27 states have opposite sign. A similar phenomenon is observed in the 
dependence of the fe-string tensions on Nc in Yang Mills theories where indeed the 1/Nc correction 
is also related with the value of the second Casimir - the so called Casimir scaling [^ . 

Finally we note that the functions are not required to vanish at x = y. Thus, unlike 

in the dipole model the solutions to the BFKL equation for ^l{x,y) are not limited to the Moebius 
invariant ones but include a larger set of solutions which are nonvanishing for the configuration 
where the transverse coordinates of two s - channel gluons coincide. This will be important for the 
discussion of Section 5. 


4.3.2 Odderon and friends. 


In the C-parity odd sector Eq. ( 4.34 ) and Eq. ( 4.35 ) lead to three decoupled equations for the 


27 - 






















functions and . For the component we have 

/ \^^u,x,y^ “ 1 “ ^u,y,x^ ( 2 /)^)] I \^^v,v,y^ (®) “ 1 “ ^v,v,x^ 


- 2 


J [Kx,y,z Kx^x,z]*^ i^^y) — ]\[ ^ ( 2 /)^) 


(4.43) 


This is the BFKL equation for the Odderon. The solutions are antisymmetric BFKL eigenfunctions. 
This solution is known as the BLV Odderon and was originally found in |M|. This Odderon is a 
’’bound state” of the reggeized gluon and the d-Reggeon. 

For the signature even 4>'*' component we obtain the very same equation but with the eigenvalues 
twice as big 

/ \^^u,x,y ^ T ^u,y,x^ (22)^)] T / \^^v,v,y*^ T ^v,v,x^ (^)22)] 

Ju Jv 

- 2 ^ [Kx,y,, - Kx,x,z] y) = ^ ^^{y, x) (4.44) 


In the t-channel the color tensor structure of this eigenstate is that of Its interpretation 

therefore is that of the bound state of the (degenerate) 10 and 10 Reggeons which is antisymmetric 
both in color and coordinate spaces. The eigenvalues for eigenstates on this Odderon trajectory 
have twice the value of the corresponding BLV eigenstates. Since for the BLV Odderon a; < 0, the 
exchange due to this second Odderon is always suppressed. As we will show in the next subsection 
its coupling in the perturbative regime is also suppressed by relative to the BLV Odderon, since 
this solution corresponds to the t-channel state containing at least six gluons. 

Surprisingly we find that the P even Odderon solution, obeys the ordinary BFKL equation: 
\2 Kx^y^z ^ S {^iV) ~ KxjZjy'ifg ~ Kz^y^x'^ s iVi^) 



+ 2 6{x-y) / Kz^u,x'^s 
Ju 


(4.45) 


- I [ 2 Kx,x,z'^s {x,y) - Kz,z,y'^s {x,z) - Kz,z,x ^ s {y, z)] = (x, y) 


N 


Since wave function is symmetric under the interchange of the coordinates the eigenvalues correspond¬ 
ing to this solution lie on the very same trajectory as the BFKL Pomeron . This Odderon 

grows with energy and dominates over all previously known C odd exchanges. We will comment on 
its couplings in the next subsection. 


4.4 Scattering of a gluonic dipole 

In this subsection we comment on the relation of the eigenstates discussed above to the cross section 
of the scattering of an s - channel gluonic dipole. The basic relation is very similar to that of the 
dipole model. Consider the evolution of the scattering amplitude of a single gluonic dipole. The 
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projectile wavefunction at initial rapidity eq. (|2.18|) consists of two gluons at transverse coordinates 
X, y in a color singlet state. The projectile averaged s-matrix S eq.(2.20|) has the form 


i 

The latter form allows us to express this wave function in terms of the two particle eigenstates of 
Hq. We use the completeness relation Eq. (|3.16|) to write: 




1 ^^ + ^ 


vr^ \x — y\^ 

+ 5*^ [tri?(x) + trii(y) — {N'^ — 1)] 




(4.47) 


Eq. ( 4.47 ) has the structure of Eq. ( 2.24 ). Let us define 7 g as 


Jq{x,y) = 


1 




— 1 TT^ |x — y|^ 


E^’’^*{x-p,y-p) 


The evolved amplitude of the scattering of the dipole off a target specified by a wave function 
at rapidity Y reads 


DSW^^iS] G+*[5] 


I q 


+ DS [tr5(x) + trS(y) - {N^ - 1)] 


(4.48) 


where the sum over q denotes summation over n and integrations over v and p. Just like in the dipole 
model this expression violates unitarity. Even if we arrange Sy(,[S'] so that all coefficients satisfy 


eq.( 2.2C ) at the initial rapidity Yq, this property is violated at higher rapidities. Some ’’probabilities” 
in eq.( 4.48| ) grow beyond unity while some become negative. The overall norm of the state remains 
normalized, Ey[l] = 1. 

The situation is in fact even more complicated, since eq.( [1.47 ) notwithstanding, gluonic dipole 
couples also to higher multi Reggeon states. Thus for example the gluonic dipole will also couple to 
’’Pomeron” containing three reggeized gluons. We discuss this in detail in the next section. These 
problems are endemic to any approximation that treats the interaction term in the RET Hamiltonian 
perturbatively and can only be solved beyond this perturbative expansion. 

Setting these problems aside for the moment we would like to ask a very specific question about 
the minimal t-channel gluon content of the Pomerons we have found in this section. To this end, we 
expand the expression for S in perturbative series. In our approach, the coupling constant enters 
explicitly only in the phase of the matrix S: 


1 


5“T'(x) = 1 + iT^^a^{x) - -{T^T^)^^a%x)a\x) + .. 


(4.49) 


~ 29 ~ 






















In perturbation theory (for small targets) the magnitude of the field a is proportional to Thus to 
find the leading perturbative content of each Pomeron we have to expand the appropriate eigenstate 
G"*" [S'] to the lowest order in a. To order only the BFKL Pomeron eigenstate contributes 

P%7S“'^(x)S^^(y) = Na^{x)a-{y). (4.50) 

To this order the scattering amplitude (we omit the terms which do not evolve with energy) is; 

‘5«2(T) = ^ 7 g(x,y) f {a'^{u) a°‘{v))T - p,v - p) (4.51) 

q Ju,v 

where the target average is defined as 

{0)t = I DSW^^[S]0[S] 


In order we get contributions also from all symmetric representations 


l\j2 \ 2 


i = 8s, 27, ^Z^ 


(4.52) 


Those come from expanding each factor of S to second order in a. Note that there is also 0{ag) 
correction to the nonlinear relation between a{x) and p{x) according to eq.(2.2). To this order this 
contributes only to the reggeized gluon exchange and therefore only to the BFKL Pomeron. The 
0{ag) contribution to the S-matrix of the gluonic dipole is 


A? 


s„.(r) = iv^Zf 


(V-Yo) 


Y^ 7 q(x,y) 


X P 


i ab 
cd 


f {a°'{u) a^{u) a^{v) a'^{v))T {u — p,v — p) (4.53) 

J U,V 


The dominant exponent at the four t - channel gluon level comes from i = 27 unless the target has a 
very particular structure which suppresses its contribution through the S'-integral in eq.(4.52). Note 
that for a projectile containing four gluons in the s-channel at order 0{ag) there is also a contribution 
due to the double Pomeron exchange. Compared to the double Pomeron exchange, the residue of 
the [27, 27] exchange is suppressed by l/N'^. 


It is also interesting that the [10 + 10,10 + I(T] Pomeron does not couple to the gluon dipole at 
0{a^). This state contains at least six t channel gluons and has a further suppression in its coupling 
strength to dilute targets. 


Let us now consider a projectile which couples to a negative charge conjugation exchange such 
as the Odderon. Note that a charge conjugation odd exchange does not contribute directly to the 
forward scattering of any projectile which is an eigenstate of C-parity. The physical states in pure 
Yang-Mills theory are indeed eigenstates of the C-parity since they are color singlets due to the 
confinement of color. In perturbation theory however the color is not confined and one is in principle 
free to consider projectiles which are not color singlets. Such projectile wave functions also do not 
have to be G-parity eigenstates. In particular one can consider a gluonic octet dipole which is 
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a superposition of the 8^ and 8s representations^. The Odderon does contribute to the forward 
scattering amplitude of such a state. Another example of such a state is the quark dipole. Although 
it is a color singlet, the dipole with a fixed transverse position of a quark and an antiquark is not an 
eigenstate of C-parity, and the Odderon was found to contribute to its forward scattering amplitude 
in 155], 156]. ^From our perspective considering such states is merely a convenient trick to make explicit 
what is the t-channel gluon content of Odderon states. 


Consider therefore the ’’state” 


S[i?] = Z\fR^^{x)Rf^Hy) (4.54) 

The lowest order perturbative expansion starts with the term of order where two t channel gluons 
couple to one leg of the dipole and the third one to the second leg. We find 

Z\f S^^ix) ^ (a“(x) a\x) a-{y) - a^y) a\y) a'=(x)] (4.55) 

The S-matrix to this order is 

5„3(y) = f {a‘^{u)a\u)a^{v)-a‘^{v)a\v)a%u))T 

(4.56) 

Where are antisymmetric solutions of the BFKL equation, see ]^^. Thus the Z^ exchange indeed 
behaves as the usual (BLV) Odderon with intercept 1 + = 0). 

For the Z~ state we can simply refer to the derivation in the Pomeron sector, since as mentioned 
above it has the color structure of —in the t-channel. Therefore just like the (10 + 10,10 + 10] 
Pomeron it contains at least six t-channel gluons. We have checked explicitly that the coupling of 
this state at the six gluon level does not vanish. 

Finally, the P even Odderon solution which arose in our calculation couples to the ’’state” 

E[P] = P“^(x) R^^iy) = i?“^(x) P^^(y) (4.57) 

Expanding to the lowest order in a we find 

S“7x) S^^iy) = — [a'^ix) a\x) a%y) + a“(y) a\y) a%x)] (4.58) 

The S-matrix to this gluon order is 

Q 

X [ {a‘^{u) a\u) a"(u) + a“(u) a\v) a\u))T E^’%u - p,v - p) (4.59) 


^Note that within perturbation theory it is perfectly legitimate to consider projectiles which are not color singlets, 
even though the cross section of scattering of two nonsinglet objects is infrared divergent. The divergence appears as 
the result of the integration over the impact parameter, while the scattering probability at fixed imp! act parameter, 
which is the focus of our discussion, is finite for colored as well as colorless objects. 
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Thus the minimal content of this Odderon solution is also three gluons and it’s coupling is not 
suppressed by powers of l/N'^ relative to the BLV Odderon. It may seem odd that this solution as 
opposed to the BLV Odderon has not been discussed in the framework of the scattering matrix of a 
fundamental dipole. The reason is that the fundamental dipole does not couple to this solution at 
any order. The fundamental dipole state is 

DF{x,y) = ^ti[E}p{x)RF{y)\ (4.60) 

The C-parity transforms it into 

C^DF{x,y)C = ^tT[Rl,{y)RF{x)] (4.61) 

Thus in order to couple to this state a C-odd exchange must also be P-odd. This excludes the 
G~^ Odderon. This is also the reason why this Odderon does not appear among the dipole 

model eigenstates discussed in Section 3 even though it is formally not suppressed by 1/N. It would 
be interesting to perform a similar analysis for more general physical states in order to understand 
whether there are any physical states that can couple to this Odderon. 


4.5 The perturbation. 


The first unitarity corrections in perturbative RFT originate from applying the perturbation Hi to 
a generic t - channel exchange in the two particle sector described by the wavefunction 


\gg) = / RfR2^nUu,v)\0) 


(4.62) 


For color singlet t-channel exchanges the state \gg) can be decomposed into components along the 
Pomeron and Odderon directions. More generally an arbitrary t channel exchange \gg) does not have 
to be a color singlet and therefore has projections onto non singlet t-channel exchanges which we 
have not studied in this paper. Quite generally however, the perturbation Hi when acting on the 
state in a two gluon sector \gg) creates a state in the four gluon sector \gggg): 


\ 999 g) = Hi \99) = Kx,y,z T“ {R^ - R^) 


afi 


Tjam hbm / r? r? \ rpb 

A\,.y ■Lh.y \^A\,y 


- 7(5 


y) |0) 


+ H 


x,y,z 


(P, - P7 


- 7(5 


^am nbm 
^z ^z 


{Ry - Rz)T^ ^lgg{y,x)\Q) 


- aS 


K^,x,zRfR 




Ky,y,zRfRz 


(P^ - Rz) 

T'^{Ry-Rz)T^ 


- ■yS 


^ggiUy x) |0) 


'x,y,z 


- 0/3 


^ggiUi x) | 0 ) 


(4.63) 


Eq. ( 4.63| ) is the most general form of the transition vertex from two Reggeons into four. To fully 
analyze the vertex in terms of the four Reggeon states we would need first to find spectrum of Hq 
in the four particle sector. This analysis is much more complicated and must involve solution of 
the BKP equation in the four gluon sector. It is beyond the scope of the present paper. We expect 
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however that the proper analysis of the transition vertex should contain, among other things the 
following elements. The ’’matrix element” of Hj between the eigenstate \gg) projected onto a single 
Pomeron (bound state of two reggeized gluons discussed above) and an eigenstate \gggg) projected 
onto two Pomerons, should reduce to the triple Pomeron vertex of Refs. [11, ^]. A generic singlet 
eigenstate \gggg) is expected to have projections onto two Odderons as well as two Pomerons. Thus 
the vertex P OO of Ref. is also expected to appear in this calculation. If the initial state \gg) 
is projected onto Odderon, we expect to find the vertex O —> ORlMj. This study is postponed until 
later date. 


5. Partonic approximation versus BFKL/BKP 


Only a small number of the multiple ” Reggeons”, ” Pomerons” and ” Odderons” that we have found 
in the previous section appear in the standard BFKL/BKP approach. It is therefore important 
to understand the relation of our approximation to BFKL and what additional physics, if any it 
contains. 


To make this connection, first let us formulate the BFKL approximation in terms of the approx¬ 
imation to the KLWMIJ Hamiltonian. In the BFKL approximation every gluon which is emitted in 
the process of the evolution is assumed to scatter on the target only once, and the scattering proba¬ 
bility is assumed to be a small parameter. This corresponds to approximating the unitary matrix R 
and its functional derivative by 


+ (5,1) 

Jikx) ^ Tf j dx- 


Denoting j dx p°‘{x,x ) by /3“, and dropping terms of order 
Hamiltonian reduces to 


^5p^(x) -I 


in eq.(2.14), the KLWMIJ 


j^BFKL 




■ (5 <5 ■ 


'6 6 ' 

.Sp'i 


SPy ^Pi. 


(5.2) 


Note that the approximation eq.(5.1) by itself as an approximation to does not 

have an expansion parameter. It becomes an expansion in only once we agree to consider the 
scattering of the projectile whose evolution we describe on dilute targets with colour fields of order 
ax = 0{as). As discussed above, after averaging over the projectile wave function, factors of 6 /6p 
turn into factors of the target field a, which has to be averaged over the target wave function. For 
dilute targets therefore, every additional factor of S/6p brings an extra power of Og- 


We can now search for the eigenstates and the eigenvalues of this functional BFKL Hamiltonian 
The calculation follows the same line as in the previous section except it is much simpler. The 


®We emphasize the word “functional” to distinguish this Hamiltonian from the commonly used BFKL Hamiltonian 
which acts on space of functions of two transverse coordinates. The Hamiltonian eq.(^.2|) generates the standard BFKL 
equation for the Pomeron in the two gluon sector, but also the whole hierarchy of the BKP equations for the higher 
”Fock space” states - see below. 
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eigenstates have the general form 


^BFKL _ T:>ab...c 

^n,X,q — 


i 


'^n,\,q{xi,X2,...,Xn) 


Xl...,Xn 


A J_ 

Xl Px2 


P\ 


| 0 ) 


(5.3) 


where projects the product of n adjoint representations onto the representation A and the 

functions ^rip,q{x\,X 2 -, ■■■,Xn) Satisfying the n-gluon BKP equations with eigenvalues ujn,\,q The 
resulting eigenstates are those and only those which are present in the standard BFKL/BKP approx¬ 
imation. In particular the only state in the one gluon sector (n = 1) is the reggeized gluon. In the 
two gluon sector the only C-even singlet exchange is the BFKL Pomeron, and the only C-odd one 
is the BLV Odderon. The nonsinglet states again are only those that can be constructed from two 
reggeized gluons. There is no trace of the ’’Reggeon” states corresponding to other color representa¬ 
tions we have discussed in the previous section, as well as of the additional Pomeron and Odderon 
states constructed from these Reggeons. 

The next question one should address is, given the eigenfunctions of the functional BFKL Hamil¬ 
tonian (BKP eigenstates) how do we calculate the scattering amplitude of a physical state ‘h. To do 
this, in principle one has to expand this state in the basis of eigenstates and evolve each eigenstate 
with its own exponential factor 


7n,A,g 

n, A, q 


riBFKL 
^n, A, q 


(5.4) 


with 

7n,A,, = {GZ%'^\^y=o) (5.5) 

This however presupposes that the eigenstates of the BFKL Hamiltonian form a complete basis of 
normalizable states. This assumption is problematic to say the least, as the states are not 

normalizable. The reason for this is easy to understand. The expansion of R to first order in Taylor 
series in eq.( |5.lD is equivalent to Taylor expansion of the eigenfunctions to the lowest order in 6 /. 
Thus the eigenfunctions in eq.( |5.3| ) being homogeneous polynomials are clearly not normalizable. 

One possible strategy to avoid this problem would be to expand the wavefunction of a physical 
state at initial rapidity in Taylor series in 6 /6p and algebraically match the coefficients to the BKP 
solutions. For example for a single s - channel gluon state we have 




^ .ah ^ ^ 

5p^{x) 2^ ^ ^ 5p^{x)6pf{x) 


+ 


(5.6) 


Then the coefficient of each power of ^ is expanded in the BKP eigenfunctions ^n,\,q- For example 
the coefficient of the linear term is given in term of the eigenfunctions of the reggeized gluon, which 


®In the rest of this section all multi-gluon BKP states will be referred to as the eigenstates of the functional BFKL 
Hamiltonian. 

^°We do not mean to imply that this is expansion of exact eigenfunctions which correspond to exact eigenvalues of 
jjKLWMij particular the spectrum of is certainly not the same as that of }jKlwmij^ latter operator 

is positive definite, while some of the BFKL eigenvalues are negative. Thus the effect of approximation eq.([S.l|) is not 
only to expand the eigenfunctions, but also to distort the spectrum of 

discussed in when the operator R acts on a symmetric function of p°'{x,x~), the path ordering in its 


definition is unimportant. We use this fact in writing eq.(5.8). 


~ 34 ~ 






are simply plane waves, 


5p°-{x) 


/ 


= Gi,a, 


q ■ 


(5.7) 


The coefficient of the quadratic term can be expanded in the BFKL eigenstates of eq.(3.13) analo¬ 


gously to eq. ( 4.47 ). Clearly only Moebius noninvariant solutions of the BFKL equation contribute 
to this expansion. The third order terms in 5/5p can be ’’projected” onto three gluon BKP states 
and so on. 

Similarly to the section 4.2, we can use in order to project the state eq.(^) onto a definite 
color representation in the f-channel . Only A = contributes to the expansion of the 

linear term while all the other representations contribute to the quadratic or higher terms only. The 
quadratic term has the structure of a ’’bound state” of two reggeized gluons at the same point in 
transverse plane projected onto a given representation A. Note that the BFKL equation for the 
channels A = 1 and A = 27 {TZj) leads to a cut and not to a simple pole in the Mellin representation. 
This is different from the contribution of the ‘exotic” colored Reggeons found in section 4.2 which 
lead to a pole structure. We will discuss the nature of the difference between these contributions 
later. 


Analogously to Eq. (5.6), for the two s - channel gluon state we would write 




+ Ap/(x 


5p^{x) 


5p^{y) 


(5.8) 


_ ^ xab rpab ^ rpcd ^ , \ |r)\ 

dp^(y)6pf(y) 6p^(x) ^ Spf(y) 


Identifying each term in the expansion with a linear combination of BKP states one obtains an 
expansion of the form 

R“'’(x) R'^"(2/)|0) = 7n, A,, (5.9) 

n,A,q 

whose evolution in rapidity is then given by eq. (M)- Note that even though some of the individual 
eigenvalues of the BKP hierarchy are positive, the scattering amplitude calculated in this way is 
expected to be unitary due to eikonalization of the multigluon exchanges inherent in the exponential 
form of R. 


Strictly speaking keeping higher order terms in the expansions eq.(5(^,5.8) is beyond the accuracy 


of the approximation eq.(5.1). Within this approximation only the lowest term in the expansion of the 
projectile wave function can be kept. More precisely one should keep the term which becomes leading 
after averaging over the target wave function. Thus for example, if the target is assumed to have only 
color singlet correlators (as any physical hadron should), the leading term is quadratic in 5/5p for 
any projectile wave function. This term describes single BFKL Pomeron exchange. Contributions 
of all higher BKP states after averaging over the wave function of a dilute target are suppressed by 
powers of ag. 

Even if one decides to keep the higher order terms it is important to realize that the procedure 
sketched in eqs.(5(^,5.9) is by no means unique. The reason is that the BKP states (|5.3D are only 


determined to lowest order in Taylor expansion. Thus for example adding higher order terms to the 


eigenfunctions is no less accurate than keeping higher order terms in the expansion eq.(5.8). On the 
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positive side, such a modification can make eigenfunctions normalizable. In particular consider the 
following natural modification: in all the eigenstates of eg. (^.31) make the substitution 


S 

(5/9“ 


-lTr(T-E) 


(5.10) 


This modification does not change the leading order terms, but now since the matrix E is unitary, 
the modified states are normalizable. With these modified states one could go back to the original 
prescription and calculate the evolution according to eqs.(^j^). This unfortunately is not quite 
kosher either, since although the modified states are normalizable, they are not necessarily orthogonal 
to each other even if they correspond to different eigenvalues. In particular a state has a 

nonvanishing overlap with a state if the wave function ..., Xn,+i) does not vanish 

when two of the coordinates coincide, Xi = xj and the overall color representations of the two states 
are the same. This is again clearly the result of the ambiguity of the higher order terms in 5/5p 
which is not parametrically suppressed within the KLWMIJ framework proper. 

This discussion is of course not original. We presented it to demonstrate that just like in the 
standard Feynman diagram approach, in the framework of the BFKL truncation of the KLWMIJ 
Hamiltonian only the one Pomeron exchange constitutes the proper weak coupling expansion, but 
also to stress that the powers of the coupling constant enter only through the procedure of averaging 
over the target. As long as we do not commit to scatter the projectile on a dilute target, that is within 
the KLWMIJ Hamiltonian proper there is no parameter which makes the BFKL truncation leading 
in any sense. The one Pomeron exchange approximation of course violates unitarity. The unitarity 
can be restored by summing contributions of higher BKP states, but this summation is ambiguous. 
In order to properly unitarize BFKL approximation truncation of the scattering matrices E in the 
Hamiltonian itself is not allowed. 


The truncation eq.(^) of the KLWMIJ Hamiltonian has two consequences, which are physically 
quite distinct. The first one is that the number of gluons in the wave function does not change 
throughout the evolution. This is obvious since each factor of 6 /6p is associated with one gluon, and 
the BFKL Hamiltonian is homogeneous in 6 /6p. As a result the eigenstates eq.(5.3) have a fixed 
number of gluons. Second, each gluon in the wave function is only allowed to scatter on the target 
once. Thus the single gluon scattering matrix 1 — i? in this approximation is always in the adjoint 
representation in the t - channel, since only single gluon exchange between any projectile gluon and 
the target is allowed. 


The complete treatment of KLWMIJ Hamiltonian will therefore bring two types of corrections 
to the BFKL approximation. The first one is due to the possibility of the gluons of the projectile 
to multiply scatter on the target. This in particular allows processes whereby a single gluon of the 
projectile exchanges color representations other than just the adjoint with the target. The second 
type of correction involves what one can call the generalized triple Pomeron vertex. It takes properly 
into account the fact that the number of gluons in the projectile wave function does not stay fixed 
through the evolution. 


/,From this perspective the partonic approximation discussed in the previous section goes beyond 
the BFKL approximation in that it allows multiple scattering of the projectile gluons. It can be 
checked directly that if we restrict the matrix E = 1 — E within the Hamiltonian eqs.(0|4.4) to have 
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only adjoint component: 


(5.11) 


R' 


ab 


-lab 


5p^ 


the partonic approximation Hamiltonian eq.(4.2) reduces directly to the BFKL Hamiltonian eq.(5.2). 
Thus clearly the partonic approximation contains all the BKP eigenstates with the same eigenvalues 
as the standard calculation. Their interpretation is just like in the BFKL approximation - ’’bound 
states” of reggeized gluons. The only difference is that now the amplitude ^ contains in principle not 
only single gluon exchange but also multiple exchanges with net color in the adjoint representation. 
The other Reggeon states in the parton approximation all correspond to t-channel exchanges in 
higher color representation and their ’’bound states” which all can occur only if each projectile gluon 
is allowed to multiply scatter on the target. As explained above these exchanges are not accounted 
for in the BFKL approximation, but instead are additional contributions which still preserve the 
total number of s - channel gluons in the projectile. 

To put things into perspective, the accuracy of including the extra Reggeon contributions into 
the calculation of the scattering amplitude is no better or worse than that of the higher BKP states. 
To begin with we have treated the variable R in the previous section as a linear variable rather than 
a unitary one, and so the states in eqs.(4.27,4.28) are strictly speaking not normalizable. However in 
complete analogy with the treatment of the higher BKP states in eq.(5.1C) we can treat R as unitary 
when calculating the norm of the eigenstates. This procedure then makes the multi-Reggeon states 
normalizable, but non-orthogonal. The non-orthogonality is of very same nature as for the modified 
BFKL/BKP states discussed above. If a wavefunction 'k(xi, ...,Xn+m) does not vanish when any m 
transverse coordinates in a state Gm+n coincide, then this state is not orthogonal to a state Gn as 
long as the two states have the same color quantum numbers. In this way the states which formally 
have n powers of R contribute to scattering amplitude of any state with m < n gluons. For example 
the Moebius noninvariant Pomeron state contributes to the scattering amplitude of a single s-channel 
gluon via a singlet t-channel exchange. The t - channel exchange amplitude in the representation 27 
will have a pole contribution from the 27 Reggeon state as well as a cut contribution from the two 
reggeized gluon state with total color representation 27 as well as from higher multi Reggeon states. 


The question whether the extra Reggeons we discussed here and their ’’composites” give impor¬ 
tant contributions to scattering amplitudes can only be decided with the knowledge of the structure 
of the target. As we have noted above, the partonic approximation is identical to the BFKL approx¬ 
imation if the target is dilute so that only single scattering is allowed. However one can imagine a 
situation in which the target is not so dilute that multiple scatterings are important but still not 
dense so that the scattering amplitude of a single gluon is significantly smaller than unity. The 
smallness in this case will not be parametric in ag, but for example simply numerical. In that case 
the other Reggeons are a priory as important as the reggeized gluon and their ’’composites” are no 
less important than the BKP states. 


To summarize this discussion, the Reggeon states discussed in the previous section are comple¬ 
mentary to the higher BKP states. Both take into account multiple scattering corrections, but the 
corrections are of different kinds. The BKP states encode corrections due to scattering of multiple 
gluons of the projectile where each gluons scatters on the target only once. The extra Reggeons 
and their composites on the other hand encode unitarization corrections due to processes where each 
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projectile gluon scatters on the target many times. In order to treat unitarization consistently in 
addition to those two types of corrections one also has to take into account the non-conservation of 
the s - channel gluon number due to the interaction term in the KLWMIJ Hamiltonian. We note 
an interesting point that some of the additional Reggeon composites we have discussed above lead 
to faster growth of the scattering amplitude with energy than the corresponding BKP states that 
contribute nominally in the same order in Og- For example as discussed above the [27,27] exchange 
grows significantly faster than the known four gluon BKP states (its overall amplitude is however sup¬ 
pressed by the factor 1/A^^). This may be an indication that the extra Reggeons and their composites 
are no less important for unitarization than the standard BKP states. 


6. Discussion 

Our goal in this paper was to relate the language of the KLWMIJ/JIMWLK evolution to that of the 
Reggeon held theory. We have shown that the KLWMIJ (or equivalently JIMWLK) Hamiltonian 
should be understood directly as the second quantized Hamiltonian of the Reggeon held theory. We 
have described precise correspondence between many elements of the standard QCD Reggeon theory 
approach and the calculation of QCD scattering amplitudes using JIMWLK evolution. We note that 
our discussion is not in terms of cut amplitudes - the objects of standard use in much of Reggeon 
held theory literature. Filling this missing link is an interesting question which should be addressed 
in future work. 

The quantum held degree of freedom in the Reggeon held theory discussed here is the SU{N) 
matrix valued held which directly corresponds to the eikonal scattering matrix of a single gluon. The 
fact that the basic degree of freedom is a local held in the transverse space makes the theory similar 
to the original pre-QCD Reggeon held theory of Gribov analyzed in |2| (and also recently in fS^ ). 
In this respect it is simpler than the Pomeron Hamiltonian considered in Section 3 (see also |2^) 
based on the dipole model, whose basic degree of freedom is the bilocal dipole held. It is however 
much more complicated than the Gribov theory in two aspects. First, the basic held is a matrix, 
and thus the number of degrees of freedom is signihcantly greater. And second, the interaction is 
nonlocal but rather decreases in the transverse plane only as I which leads to the violation of the 
Froisart bound Both these features are related to the specihc nature of the QGD interaction. 
The proliferation of degrees of freedom is simply due to the fact that the basic partons of any QGD 
state are gluons, and the nonlocality of the interaction is the direct consequence of the perturbative 
masslessness of the gluon. One would hope that proper account of conhnement can lead to an effective 
theory, valid at large transverse distances jxj > 1/Aqcd which operates directly only with scattering 
amplitudes of gauge invariant states and is local on these transverse distance scales. Such a theory 
should follow from the fundamental QGD RFT by operator product expansion. It is however unclear 
to us whether the number of degrees of freedom in such a theory should remain finite. At any rate 
this question is far beyond the scope of the present research. 

In this work we have considered perturbative approach in the framework of the QGD RFT and 
have found several eigenstates of the free Hamiltonian. In particular we have found all eigenstates 
in the one particle channel and color singlet states in the two particle states. 
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In the one particle sector we find the reggeized gluon as well as a set of additional Reggeons 
corresponding to the exchanges in the color representations 85 , 10, 10, 27 and TZj. In the color 
singlet two particle sector we find states which are naturally interpreted as bound states of these 
Reggeons. In particular we find the standard BFKL Pomeron which is the ’’bound state” of two 
reggeized gluons, and also Pomerons that correspond to symmetric bound states of two 85 , two 27, 
two 10 + 10 and two TZ^ Reggeons. All these correspond to charge conjugation even exchanges in the 
f-channel and have positive signature. 

An interesting observation here is that the [27, 27] Pomeron exchange grows faster with energy 
than two BFKL Pomerons. Also the large Nc correction to its intercept is 0{1/Nc) and not 0{1/N^) 
as could be naively expected. The reason for this large correction is that in the Aj, — 00 limit the 
[27,27] and [ 77 . 7 , 777 ] Pomerons are degenerate and so the 1/iVc expansion is actually a degenerate 
perturbation theory. 

We also analyzed the C parity odd sector. Here we find three sets of eigenstates. One is the 
signature odd bound state of the reggeized gluon and the 85 Reggeon, which has all the properties of 
the BLV Odderon. Another set of eigenstates can be thought of as an antisymmetric bound state of 
10 and 10 Reggeons. These eigenstates are charge conjugation odd and parity odd. They are quite 
distinct from the standard Odderon however in that they have even signature. The intercept of this 
trajectory is unity and so these exchanges lead to cross section that does not rise with energy. Finally 
the third set of states in the Odderon sector is somewhat mysterious. The eigenstates of this set have 
negative signature but positive parity. Their wave functions are symmetric under the exchange of 
the transverse coordinates of constituents (the reggeized gluon and the 8 s Reggeon). The intercept 
of this trajectory is equal to that of the BFKL Pomeron. 


We have analyzed the t - channel gluon content of these eigenstates and have found that the BFKL 
Pomeron (unsurprisingly) contains at least two gluons, the [ 85 , 85 ] [27, 27] and [R 7 , Rj] states contain 
at least four gluons, the two signature odd Odderons contains at least 3 gluons while both [10,10] 
bound states have at least 6 gluons. Thus perturbatively the couplings of [ 85 , 85 ], [27,27], [ 777 , 777 ] 
exchanges are suppressed by and of the 10 exchanges by af relative to the BFKL Pomeron. 


7 ,From the s-channel scattering point of view this means that these higher Pomeron exchanges 
are due to multiple scatterings. Thus for example in order to exchange the [27, 27] Pomeron each 
projectile gluon has to scatter at least twice from the target as is clear from the derivation in eq.( 4.52| ). 
7 ,From the point of view of the target this means that the target field must be dense enough so that 
at a given point in the transverse plane there is significant probability to find a multigluon state in 
the representation 27 of the color group. Such scattering is suppressed if the target is dilute, that is 
if the probability to find an extra gluon at a given point in transverse plane is 0{as). However in 
the high density regime this suppression disappears, since the probabilities to find an octet or a 27 
are equal when the average color charge density is large. 


We have analyzed the relation of the spectrum we have found in partonic approximation to 
KLWMIJ with the standard BFKL/BKP approach. We have shown that all the BKP states do 
appear as eigenstates on the partonic Hamiltonian. The extra states we found are outside the 
BFKL/BKP framework, as multiple scatterings of a single projectile gluon are a fortifiory excluded 
from the BKP approximation. These are therefore physically distinct unitarization corrections. Since 
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they are related to multiple scatterings of a single projectile gluon, it is natural to expect that they 
are suppressed by powers of 1/A^ relative to the corresponding BKP states. This is indeed what we 
find. For example the contribution of [27, 27] Pomeron which has an intercept significantly greater 
than any of the four gluon BKP states, has a 1/A^^ suppression factor relative to the leading BKP 
state. 


Next we discuss some issues related to the structure of the symmetries of the RFT. As we have 
already remarked, the theory has a SUl{N) 0 SUji{N) symmetry as well as the discrete charge 
conjugation symmetry R CRC and the signature symmetry R —> R^. Expansion around R = 1 
breaks spontaneously SUl{N) 0 SUl{N) down to SUv{N) while both discrete symmetries remain 
unbroken. The breaking of the symmetry leads to the appearance of the Goldstone boson, which 
is the reggeized gluon. Since the intercept of the other Reggeons is also unity it suggests that the 
Hamiltonian may have higher symmetry. The most general transformation that one can define on 
the space of adjoint unitary matrices is GL{N‘^ — 1,R). This is obvious since R can be written in 
terms of — 1 real parameters as 7?“^ = [exp{^'^r'^}]“^ and the parameters can be transformed by 
an arbitrary real linear transformation. The SUl{N) 0 SUr^N) transformation is a small subgroup 
of this group. It would be interesting to check whether or not the KLWMIJ Hamiltonian is invariant 
under a larger subgroup of GL{N‘^ — 1, R) which is also broken by 72 = 1 configuration. There is also 
a possibility that the appearance of many massless particles has a similar origin as in supersymmetric 
theories where additional symmetry appears if one artificially extends the Hilbert space of the 
theory. We note that it has been previously suggested in Ref. [^] that the degeneracy between /- 
and d- trajectories may be related to a variant of supersymmetry [pT|]. 


Another issue related to symmetries is the status of the conformal symmetry in the framework 
of the second quantized RFT. It is well known that the BFKL equation is conformally invariant on 
the Moebius invariant sector, that is on the set of functions F{x, y) which vanish when x = y\82, 0 ]- 
The nonlinear Kovchegov equation j^j is also conformally invariant. In fact it is easy to see that the 
full KLWMIJ Hamiltonian is conformally invariant if we were to substitute the Weiszacker-Williams 
kernel in eq.(^) by the dipole kernel 


Kx,y,z M{x,y,z) = 


as {x - yf ^ 
{z — xY{z — yY 


Q!.g 


1 

{z — x )2 


1 

{z-yY_ 


( 6 . 1 ) 


Under the inversion transformation (with x = xi+ ix 2 etc.) 


R{x) 72(l/x), 


* —2 


{x*x) 


6 R{x) 6R{l/x) 

The dipole kernel can be written as 

Mix,y,z) = {z*z)~‘^ M{l/x,l/y,l/z) 


( 6 . 2 ) 


(6.3) 


Changing the integration variables from x,y,z to l/x,l/y,l/z we see that the Hamiltonian eq.( | 2 . 1 lD 
with the dipole kernel is invariant under this transformation. The full conformal invariance follows 
from the dilatational and inversion invariance. 


The KLWMIJ Hamiltonian eq.( 2.11j ) is not invariant under inversion, as it can be easily checked 
that the property eq.(|6.3|) is not shared by the Weiszacker-Williams kernel K{x,y,z). However as 
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has been extensively discussed in the literature (see or |^) if one limits oneself to consideration 
of color singlet impact factors and color singlet exchanges, the kernels K and M are interchangeable. 
Formally speaking, taking the difference between the corresponding KLWMIJ Hamiltonians one finds 
terms proportional to either 


or 



These are the generators of the SUl{N) and SUi{{N) transformations respectively. Thus when acting 
on states which are SUl{N) (8) SUfi,{N) invariant the two Hamiltonians are identical. 


We therefore conclude that although the KLWMIJ Hamiltonian is not conformally invariant, 
part of its spectrum is. Note however that the perturbative approach described in the present paper 
breaks spontaneously the SUl{N) (g) SUr{N) symmetry. This type of perturbation theory would 
give different results for the two Hamiltonians. In particular the dipole variant of the KLWMIJ 
Hamiltonian has no virtual terms of the type eq.(4.4). Its one particle spectrum is therefore trivial 
- any one particle state is its eigenstate with zero eigenvalue. On the other hand the spectrum of 
the Pomerons and Odderons discussed above remains unaffected. This may have to do with the fact 
that even though the perturbative expansion is not SUl{N) 0 SUr{N) invariant, some eigenstates 
nevertheless do have the full symmetry of the Hamiltonian. 


We note that we have not made a serious attempt to formalize the approximation scheme em¬ 
ployed in this paper. In particular it has two important elements. First we have split the RFT 
Hamiltonian into the free and the interaction parts according to eqs. (|4.2| ) and ( [4.5D . The idea behind 
it is to separate H into a homogeneous and an inhomogeneous part. The homogeneous part preserves 
the number of s-channel gluons throughout the evolution, while the inhomogeneous term increases 
this number. 


However such a separation is not unique due to the fact that Ris a. unitary matrix. Our choice 
of the split was guided by two criteria: we required that the free and the interaction parts are 
separately ultraviolet and infrared finite. The infrared finiteness is straightforward, it simply asks 
not to emit new gluons far away from the existing gluons with probability greater than the square of 
the Weiszacker-Williams field. The origin of the requirement of the ultraviolet finiteness in physical 
terms is the following. Calculating the cross sections using the homogeneous part of the Hamiltonian 
corresponds to the partonic approximation, which treats the scattering of the projectile gluons as 
completely independent. However physically it is clear that if two gluons in the wave function are 
separated by a distance which is smaller than the correlation length of the fields in the target, 
such partonic approximation is not valid. In this case the two gluon system scatters like a single 
object in the color representation that combines the two gluons and the color correlations can not 
be neglected. In particular if this two gluon state emerges from a single gluon through a step in 
the evolution, such state should scatter as a single octet, and not as an independent product of two 
octets. For separations greater than the correlation scale the partonic approximation makes perfect 
sense. However our approach is designed for an arbitrary target and the split into Hq and Hj can 
not depend on ! the characteristic of the target. The physical alternative then is to choose such a 
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split which suppresses the emission of a gluon at a neighboring point without introducing a scale, 
that is in a dilatationally invariant way. Technically this is ensured in eqs.(4.3,^!4) by the presence of 
the factors R{x) — R{z). This is the only possibility we found to satisfy the condition of ultraviolet 
finiteness, but we have not proven that it is unique. 


The second important element of our approach was the diagonalization of Hq. To do this we 
had to differentiate functions of R. Since i? is a unitary matrix, this has to be done with care. 
The complete RFT Hamiltonian preserves the unitarity of R. This means that when acting on a 
wave function which vanishes for nonunitary R, the Hamiltonian transforms it into another wave 
function with the same property. In this situation one can substitute the constraint differentiation 
with respect to unitary R by an unconstrained one as it does not change the action of H on physical 
states. However this property is not satisfied by Hq and Hj separately. In our calculation we have 
disregarded this subtlety. A more careful approach would be to introduce the Lagrange multiplier 
that imposes the unitarity of the matrix R. Although we do not believe that this will affect our 
results, it is a very important question worth detailed study. In particular there is in principle a 
possibility that some of the states we found may be pushed into the unphysical part of the Hilbert 
space once the constraint on the R matrix is imposed. The prime suspect here could be the parity 
even Odderon, since it does not couple to simple color singlet states. 


As we have discussed in Section 4, the perturbative treatment of H violates unitarity. This leads 
to appearance of negative probabilities as well as to the unbounded growth of the cross section. In 
this context the Pomeron states are tachyonic, since they have negative eigenvalues at low momentum 
transfers. In fact the Hamiltonian Hq is not bounded from below and it is clear that such tachyons 
exist in sectors with arbitrary number of s - channel gluons with ever increasing intercepts. The 
cure for this can only come by including the interaction Hj nonperturbatively. We certainly believe 
that the JIMWLK/KLWMIJ Hamiltonian does define a unitary theory. Numerical studies of both 
the BK equation and of the full JIMWLK equation lead to this conclusion [^ ]. It would be very 
interesting to devise a nonperturbative approach to RFT, perhaps along the lines of a mean field 
approximation with {R) = 0. 


A state with {R) = 0 corresponds to the black disk limit. The SU{N) ® SU{N) symmetry is 
restored in such a state. In physical terms the restoration of the symmetry means that the color 
index of an incoming gluons is not correlated with the color index of the outgoing one. Thus the 
color is completely randomized during the interaction with the target. It is usual in the quantum 
field theory that the restoration of symmetry is associated with condensation of the excitations of 
the ordered phase. In our case those excitations are Pomerons. We can think therefore of the black 
disk limit as condensation of Pomerons similar to the old ideas of (see also [^| for the renewed 
interest in this problem). 


This analogy is not perfect. There is no phase transition in RFT, since there really are no two 
distinct phases as a function of some external parameter like temperature or coupling constant. The 
only parameter in the RFT Hamiltonian eq.( |2.11 ) is but it enters as an overall multiplicative 
constant and thus its value can not affect the phase structure. It only affects the overall scale of the 
energy and therefore the speed with which a given state evolves towards the vacuum. So one and 
the same phase must be realized at all values of ag. The Pomeron condensation therefore occurs 
always. This is consistent with the fact that as excitations they have tachyonic nature already 
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in the perturbative expansion. In principle the situation might change when the additional terms 
responsible for Pomeron loops are included (see discussion below), although we would still expect 
the black disc limit to be the ground state for all values of ag at least as long as nonperturbative 
confining effects are not accounted for. 

A similar situation is encountered in the Gribov effective Reggeon field theory |^, where the 
nonperturbative account of the Pomeron interactions leads to a gray disc limiting behavior with 
the growth of the total cross section which satisfies the Froissart bound. There are however many 
differences between the two cases. Apart from the ones already mentioned above, one very important 
difference is the fact that the KLWMIJ Hamiltonian which has been the basis of the present work 
violates t-channel unitarity as it does not include Pomeron loops. 


There has been a lot of activity recently in an attempt to incorporate the Pomeron loops into 

1 , 13 , g, [ 


the JIMWLK/KLWMIJ formalism 


In spite of some progress 


in this direction the complete evolution kernel is still unknown. We therefore still do not have a 
completely unitary RFT. One could try naively to restore the f-channel unitarity by symmetrizing 
the Hamiltonian with respect to R and This would be analogous to the Gribov theory which 
is symmetric under the interchange of the Pomeron field and its conjugate momentum. This idea 
would lead us to add terms of the type to the Hamiltonian. However these terms 

correspond to the process of recombination of three s - channel gluons into two as a result of the 
boost of the projectile wave function. The s-channel gluons however do not recombine, but rather 
are created by the evolution (see [^| for discussion). The nonlinearities in the JIMWLK evolution 
for example correspond to the decrease in the emission probability of an extra gluon, but not to 
the appearance of recombination in s - channel. In Refs. we have derived an extended 

version of the JIMWLK/KLWMIJ equations (JIMWLK+/KLWMIJ+), which account partially for 
the effects of Pomeron loops. Similar results were obtained in the path integral approach in 
|5l| . This work suggests that the t-channel unitarity corrections should occur as higher powers of 
This should have an effect of generating n ^ m {n > m) Reggeon transitions alongside the 
m ^ n transitions present in KLWMIJ. It remains unclear whether the effects accounted for in the 
JIMWLK+/KLWMIJ+ equations are indeed sufficient to restore the t-channel unitarity of the RFT. 
It would be interesting to understand in general whether the selfduality condition of Ref. |g| ensures 
the t-channel unitarity. The f-channel interpretation which emerged naturally in the present paper 
should be helpful in addressing these questions. 


Finally we note that in this paper we considered RFT in gluonic sector only. Although in the 
leading eikonal approximation only gluons are created in the evolution, it is not difficult to consider 
initial states which contain quarks. To this end essentially all one has to do is to rewrite the left and 
right color rotation generators in eq. (|2.1I ) in terms of the fundamental rather than adjoint matrices 
Rp. One can then repeat the perturbative analysis of the present paper without any additional 
complications. It may be interesting to consider the quark dipole scattering as it is more directly 
related to physical process (such as DIS) than the gluon dipole. 
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A. Appendix: Evolution of 


In this appendix we consider the evolution of the four gluon t - channel state coupled to a single dipole 
(s - channel) projectile. The discussion is in the framework of the leading large N approximation. 


The state with n gluons in t-channel is defined by eq.(2.25). We consider two such states - with two 
and four gluons. We also project the gluons pairwise on color singlet states as appropriate for the 
large N limit. We thus consider 


D^{Y){xi,X2) = 


1 

A 


Dpp^{xi)p°‘{x2)Wy[p] 


and 


D^{Y){xi,X2.,X:i,X4) = 


A2 


Ap/ (xi) / (X2 ) /o'* (X3 ) (X4 ) Wy [p] 


(A.l) 


(A.2) 


There is one subtle point in this definition which turns out to be important for the calculation 
of the initial condition. Namely when any two transverse coordinates are equal (xi = X 3 , etc.) 
the appropriate charge density operators do not commute with each other. In this case one has 
to remember that in the integral in eq.( |A. 2 | ) the charge density variable is also endowed with the 
longitudinal coordinate x~, which keeps track of the ordering of the operators p“(xi) —> p“(xi,xj~). 
This has been discussed in great detail in [^] where we have also developed a technique to calculate 
this type of averages taking the noncommutativity into account. We do not give details here but 
interested reader can consult |^]. 

The function is the one considered in |l^. We choose to concentrate on it rather than on a 
function completely symmetric with respect to all four gluons | 12 | because the color algebra here is 
marginally simpler. The evolution with rapidity is given by 


—D^ = - 
dY N 


A 

dY^ 


= 


1 


Dpp-{xr)p-{x2)x^^^^^^[p, 

[ Dpp%xi)p-ix2)p\xs)p\x,)x^^™-^[p, j^]WY[p] 


(A.3) 


. Eor a ’’virtual photon” projectile state with wave function P'^{u,v) the initial projectile weight 
functional is 

'^[p] = f P'^ {u,v)r{u,v)6[p] (A.4) 

Ju,v 

with r{,v) defined in eq. (|3.1| ).). Eor such a projectile the initial conditions for eq.(A.3) are 

1 


D^iY = 0) = f4o(^i,X2) = 


N 


P'^{u,v) / Dpp°-{xi)p°-{x2)r{u,v)5[p\ 


I 


= / P'^{xi,v)5xi,x2 - P^ixi,X2) 


(A.5) 
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= 0 ) = i:»o(xi,X 2 ,X 3 ,X 4 ) = ^ J P'^{u,v) j Dpp'^{xi)p'^{x2)p^{x‘i)p^{xi)r{u,v)5[p] 

1 


DQ{xi,X‘i)5, 


• u,v 

'XI,X2 '-'XZ,X4 + Dq(xi, X 2 ) ^Xl,X3 ^X2,X4 + Dq(xi, X 2 ) <5x1, 0:4 <5x2,X3 


(A.6) 


- DI{xi,X 2 ) ^X 2 ^XZ ^X2,X4 - Dq(xi,X 2 ) ^Xl,X 3 ^Xl,X4 


- Do(xi,X3) 

^Xl,X2 ^Xl,X4 - Dq(xi,X4) 

<5x1,X 2 <5: 


XI,X3 


It is straightforward to verify explicitly by expanding r(u, v) in powers of 6/6p that the functions 
and Dq are the standard impact factors for coupling of two and four gluons to a dipole projectile. 
The explicit relation between Dq and Dq is the same as that of m- We stress again that in the 
derivation of eq.(A.5) it is important to account correctly for noncommutativity of operators p°' at 
the same transverse coordinate [41|. 

To write explicitly the evolution equations for and we use ^^lwmij dipole limit 

eq.( |3.3|) and act with the kernel to the left in eq.(A.3)^^ Since most on fourth power 

of p, the nonvanishing contributions come only from expanding operators r {d — 1) in the interaction 
part of the kernel to order {5/5p)‘^. The leading large N contribution of the free part comes from 
the same order of expansion of r. Thus for the purpose of the present calculation the dipole creation 
and annihilation operators are 

r r n 2 


d\x,y) = 


AN l6p‘^{x) 5p°-{x) 

and the kernel can be written as 
1 


d{x,y) = 


X = 


4iV2 
2 


J x,y,z 


6p^{x) 6p^{z) J \Sp^{y) 6p^{z) J AN \6p^{x) Sp^{z) J \6p^{y) 5p^{z) 


(A.7) 

(A.8) 
2" 


Now acting with this kernel in eq.(A.2) we obtain 

d f 

—D‘^{xi,X 2)= K{xi,X2-,yi,y2)D^{yi,y2) (A.9) 

Jyim 

and 

d f 

—D'^{xi,X2,X3,X4) = / [Kixi,X2-,yi,y2)D'^iyi,y2,X3,x4) + K{x3,X4;yi,y2)D^ixi,X2,yi,y2) 

dyi,y2 

- V{xi,X2,X3,X4,yi,y2)D'^{yi,y2)] (A.IO) 


where K{xi,X 2 ',yi,y 2 ) is the nonforward BFKL kernel and V is the so called three Pomeron vertex 
V{xi,X2,X3,X4,yi,y2) M{xi,X3,X2)dy^ ,Xi<52;2,X3<5x2,X4 + M(xi,X4,X2)(5yi 

,Xl<5y2,X4<5x2,X3 

+M{x 2 ,X 3 ,Xi) 6 y, 

,X2^1;2,X3^X1,X4 + M{x 2 ,X 4 ,Xi) 6 y^ 

,X2 ^y2 ,X4 dxi,X3 (A.ll) 


^^Note that the full KLWMIJ Hamiltonian beyond the dipole limit contains additional terms in the expansion to 
fourth order in derivatives. The additional terms however do not contribute when acting on and D'* as dehned 
above, see 
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In eq.( A.ll ) we have not written explicitly terms proportional to Sxi,x 2 and/or Sx 3 ,x 4 since they vanish 
when multiplied by the Pomeron eigenfunctions of the dipole model and thus do not contribute to 
solution of eq.( A.IOD , see for discussion. 

Eq.(A.9) is clearly the BFKL equation. Solution of the nonhomogeneous eq.(A.lO) can be written 
in the form 

D^(Y) = Df(Y)+DUr) (A.12) 


Here the irreducible piece Df is defined as the solution of eq.( |A.lCi ) with initial condition Dj{Y = 


0) = 0, while the reggeized piece is dehned as the solution of the homogeneous part of the 


equation (eq.( A.l^) excluding the last vertex piece) with the initial condition D‘^{Y = 0) = Dg. The 


piece is precisely the reggeized piece in the standard terminology. Its dependence on energy is 
that of the single Pomeron exchange rather than the double Pomeron exchange. The reason is that 
the initial condition Dq of eq. (|A.5| ) is given by a linear combination of Dq due to the fact that the 
dipole projectile has only two transverse coordinates. On the other hand the homogeneous part of 
the evolution equation eq. (|A.10 ) has the property that it propagates such an initial condition with 
the BFKL kernel. 


The reason why the reggeized piece appears in the present analysis but does not appear in 
the analysis of the dipole model presented in [l^ is the following. The authors of |l^ analyzed 
the equation for the inclusive double dipole density rather than for the four gluon correlator. The 
two quantities although similar, are not completely equivalent. In particular it is obvious that the 
inclusive double dipole density vanishes in the initial single dipole state. Thus even though in the 
large N limit it satisfies the same evolution equation as the four gluon correlator, the solutions of 
the two equations differ precisely by the reggeized piece. 

We thus see that our approach reproduces the standard splitting of the gluon four point function 
into the irreducible part and the reggeized part at least in the large N limit. Comparison of this 
aspect of the two approaches in more general setting is beyond the scope of this work. 


B. Appendix: Derivation of the eigenvalne equations 

In this Appendix we give details of computation of the left hand side of eqs. ( 4.33 , 4.34 ). We start by 
acting with the Hamiltonian Hq on the functions We begin with and then repeat 

the analysis for the function \ Acting with the ’’real” part of the Hamiltonian Hq we obtain 

HRG+i++) = -2ku,v,z {t:^[Ru - R^r^T^kiRv - + [Ru - r|, [i?, - R.rT^a 

- [Ru - [Rx - R^m^ - [Ru - Rzk^ [Ru - Rzrr-,] (b.i) 

Similarly applying the “virtual” Hamiltonian Hq gives 

G+(++) = A f Ku,u,zPii^^i{u,v) {{T~^T%d[Ru - R,f - iV [i?„ - R.T^} lP\u) 

Ju.v.z ^ ^ 


m 


(B.2) 
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The next step is to differentiate Eq. (BT) and Eq. (B^) with respect to 




SRT 6K 


V 


tjR /^+ (H—h) I '~pCL rpa Tfi f3S I rjia rpa -pi bd rjia n-ia pi (Bd rjia rjia pi dS I 

0 ^ 1 aa ac ' ^ /Bb ^ Sd ay ^ aa ^ 5d ay ^ 0d ^ay ^ aa j 


^l{x,y) - K^^:,^y^l{x,z) - z) +25{x-y) / Kz^u,x^\{z.,u) 


L 




'L 


2Kx,y,z %{x,y) - Kx,z,y^\{x,z) - Kz,y,x^\{y-,z) + 26{x-y) / Kz^u,x%{z.,u) 


(B.3) 


The last equality follows by using the crossing properties of the projectors Eq. ( |C.18 )-Eq. (C.28). 
The contribution of the virtual part is 

SR-yx) SR'ryy) ° ^ / 

® / [(Py,j/,2: T Rx,x,z) ^s(Xj y) R-z,z,y ^s(Xj z) Kz,z,x ^siV^ z)\^ 

J z 


= AN 


Cj^pk _ pi 


0S 
J ay 


^ l2K,,,,,A>i(x,y) - K,,z,y^i(x,z) - K,,z,x'I'Uy, z}] (B.4) 


where again Eq. ( C.18 )-Eq. (p.2^) were used to get the last equality. 

We repeat the same algebra but now applied to the function \ The “real” and “virtual” 

terms read 

) = -2K^,,,z {Ty[Ru - RzV^Ty[R, - P,]“ + [P„ - P,]“/^r|,[P, - 
- ry [P„ - RYY^ [P„ - - ry [r, - Ry^<^ [p„ - (b.5) 


ny G^(—) = A [ K,^,,zZty^^{u,v) {(T~^T\d[R, - R,f - N [R, - R^\u) 

J u,v,z ^ ^ 


(B.6) 


Differentiating Eq. ([B.5|) and Eq. (|B.6| ) with respect to obtain 

<5 6 


5Rf dPf 


6 R°‘P{x) SWRy) 

ZjR^^i -) _ J 'T~'^ I 'T'^ 'T'^ ry-^bd rpa rpa ry-^ (Bd rpa rpa ry^dS I 

0 ^ aa ^ cy ^ ac ' (Bb ^ 8d ^ ay ^ aa ^ 6d ^ ay ^ (Bd ^ay ^ aa j 

(g) 2 y [2Ka;^y^z^'^{x,y) - K^^z^y^^{x,z) - K:z,y,x^^{z,y)\ 

= [2K^^y^,A>^{x,y) - K^,z,y^^{x,z) - K,^y^^^^{y,z)\ 

(B.7) 
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We again rely on the crossing properties of the projectors Eq. ( C.18 )-Eq. ( p.31 ). Einally the contri¬ 
bution of the virtual part is 


<5 (5 


= 4 {{T'^T%d 




±I3S 


} 


Zi j \jyt^y,y,z Kx^x,z) Kz,z,y i ^') z,z,x , y')~^ 

= _2(3ti)/2^^±^^5 ^ [2Kx,x,z^^{x,y) - Kz,z,y^^{x,z) - Kz,z,x^^{z,y)] (B.8) 

The same algebra applied to the function G~ ^ gives (we note that = Z~^^) 

zjR r^— (—h) o J r p p 17^ r p p _i_ rp p r p p 

-^0 ^ — —^I\u,v,z — ^z\ ^ckl-^v ~ ^z\ + [-^u “ ^z\ ^ l3b [-^v “ ^z\ J^nd 

- [Ru - RzV^ [K - RzV TL - ^7 [Rv - RzV'^ [Ru - Rzr T^b} ^“at v) (B.9) 


<G-(-+) = 4 [ 

J U,‘ 




{u,v) [{T-T\d[Rv - Rzf - N[Rx- RzY^^] IP\u) 

(B.IO) 


Differentiating Eq. (|B.9|) and Eq. (|B.10 ) with respect to <5 _r4 (y) obtain 


a/3 .575 


6Rx^ 5R 


y 


TZtR /^— { h) j rpa rpa ry —c5 I rpa rpa ry —7^ rpa rjia ry —7(3 rpa rjia ry—aS I 

-"0 ^ ~ ~ aa c'y ^ af} P ^ 0b ^ 5d ^ ah ~ ^ aa ^ 5d afi ~ 0d ^ a-y ad j 

Z> 2 J ^ 2 Kx^y^z^s {XiV) Px,z,y^s (x^z) Kz,y,x'^s (2-) 2/)] 

^ ^ a 0 j ^ '^a:,j/, 2 'hg (x, y) — Kx^z,y 'I'g (x, z) — Kz^y^x^s {Vi 

5{x-y) / Kz^u,x^^{z,u) 

J U 


= 2 
+ 2S(x - 


(B.ll) 


The contribution of the virtual part is 
5 <5 


diif 5 r:i^ 


G"("+) = 4 - N Z~'^^ 


a/3 


0/3 


} 


[ [iRy,y,^ P Px,x,z) s {x^y) Kz,z,y'^s (XjZ) Kz^z,x^s iz:,y)] 
J z 


= -2NZ 


— 7(5 

0.(3 


[2Kx,x,z^s 7iy) - Kz,z,y^s 7)^) - Kz,z,x^s ?/)] (^-12) 


C. Appendix: Projectors in SU{N) 

In this Appendix for convenience of reference we collected some standard formulae for the SU{N) 
group and its representations. In particular throughout our calculations we have used extensively 
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the decomposition of the product of two adjoint representations into a direct sum of irreducible 
representations. We give definitions of the relevant projectors and summarize their properties. This 


is taken almost entirely from Ref. [57| and we use notations of that paper. The same algebra has 
also been worked out in |74]. For more details we refer the reader to Refs. [57, 74, The SU{3) 

results can also be found in [H, ^]. 

The 5f7(A^)-generators in the fundamental representation, t°',a = 1... — 1 define the familiar 

/- and d-tensors through 


+ 2 (^dabc + ifabc^t'^ , 

It follows that 

ifabc = , dabc = . 

The / and d tensors obey the Jacobi identities 

'^fkamifmbl '^fkb mifmal — ^fabmifkml j 


(C.l) 


(C.2) 


(C.3) 


and summation rules 


fkamdmbl dfcfc mfmal — fabmdkml ? 


faij fbij ^ ^ 


daij dfyij — 


Ar 2 _4 


N 


faij d})ij — 0 . 


(C.4) 


(C.5) 


N 

fiajfjbkfkci "irfabc ; 

..... 


fiajdjbkdkci 

diajdjbkdkci 


Ar 2 _ 4 

2N 
_ 12 


2N 


fabc 
dabc ■ 


(C.6) 


The SU{N) generators in the adjoint representation = —ifabc satisfy standard commutation 
relations [r“,r^] =ifabcT''. 

We now consider the product of two adjoint representations of the SU{N) group and its decom¬ 
position into a direct sum of irreducible representations. 

{N^ - 1) X (1V2 - 1) = 1 + (iV2 _ 1)^ + (^2 _ 

- A){N^ - 1) ^ i{N^ -A){N‘^ -1) 


^ N^{N + 3){N -1) , N^{N-3){N + 1) 
^ 4 ^ 4 

= 1 + 8a + 8s + 10 + To + 27 + 7^7 . 


(C.7) 
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symmetric 


antisymmetric 


Name of rep. 1 85 

Dimension 1 — 1 

Casimir 6*2 [7^] 0 N 


27 

Ar2(Ar+3)(Af-l) 

4 

2(iV + l) 


777 

Af2(Ar-3)(Af+l) 

4 

2(A^- 1) 


8a 

iV2 - 1 
N 


10 

(Ar2-4)(Af2-l) 

4 

2N 


10 

(7V2-4)(7V2-1) 

4 

2N 


Table 1: Properties of Multiplets 


We follow the nomenclature of |57] in denoting the representations by their dimensions in the SU{3) 
case. The symmetric representation TT-y does not exist in SU{3). In Ref. |74] this representation is 
denoted as 0 . 


The S'C/(7V)-projectors into the singlet as well as the two adjoint multiplets have manifestly the 
same form as their well-known N = 3 counterparts: 


mfd = (C.8) 

P[^A]cd = J^fabkfkcd = -J^'ifabk'ifkdc (C.9) 

= J^/abkdkcd = [DsC (C.IO) 

The projectors onto the antisymmetric 10 and 10 are 

pmfi = 5 (-^25 - p[^a\£, - iis‘). (c.ii) 

The remaining representations are symmetric and the corresponding projectors are 
Pmfd = ^ ((^ + 2)5^ - (TV + 2)(TV _ 

- i(iV - 2)(iV + 4)P[85]“^ + + [D^O] (C.12) 

P[P7]fd = ^((^ - 2)5“] + (TV - 2)(TV + 1)P[1]“] 

+ l(iV + 2)(iV - 4)P[85]“] - y ([A]“] + • (C.13) 

The tensors A, Y, S, and D are dehned as follows 

5“] = 2 + Sadhc^ , Acd = - (dachd “ ^adh^ ■ (C.14) 


^ 


(^dadkfkbc T fadkdkbc 


— dackdkbd j [VliiJcci — j [TT^sJcd — ^afefc^fccd 


(C.15) 

(C.16) 
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Apart from the complex, but hermitian structure iYg = Ps[10] — Ps[l0], which already appeared 
in the decuplet projectors, the full set of color-singlet four-gluon states includes two more complex 
hermitian tensor structures 



The projectors satisfy the Fierz-type identities: 


(C.17) 


9 

Pt[n,] = j2ciPsm, (C.18) 

i=l 

where the t-channel projectors are defined as: 

Ptmtd^Psmz- (C. 19 ) 

The Fierz relations are therefore 

= (C.20) 

i=l 


with the crossing matrix Cf (Fig. is now obtained as 


cl = 


mj]td • 


cd 

'^\ab 


m. 


ac 

-j\bd 




'^iab 


m&mii dim[P, 

Note that eq.( C.20 ) includes all nine tensors (i.e. P and Z^). 
Below we present several identities which we found useful. 


(C.21) 


{T‘^T^)m = NXjP[]Zj] 


ac 

bd 


where Xj is defined as the second line of the matrix C: 

A, ^ 

If Di as the dimension of f-th representation we have 

Ci^ = Xi — 

It is useful to introduce the matrix C (table 


i=i 


which obeys 


(C.22) 


(C.23) 


(C.24) 


(C.25) 
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Name of rep. 

1 

8a 

8s 

10 

To 

27 

Pr 

1 

8a 

85 

10 

To 

27 

Pr 

0 

1 

0 

0 

0 

0 

0 

1 

TV^-l 

1 

4 

1 

4 

0 

0 

N + 3 

0 

1 

4 

1 

4 

1 

4 

1 

4 

0 

0 

0 

0 

1 

N‘^-A 

1 

4 

1 

4 

N + 3 

0 

0 

1 

-4 

1 

4 

1 

4 

Ar + 3 

0 

1 

w 

0 

(Af + 1) (iV-2) 
4lV^ 

(Af + 1) (Ar-2) 

4lV^ 

(7V + 1) 

0 

1 

w 

0 

(iV-i)(Ar + 2) 
4lv^ 

(Ar-l)(Ar + 2) 
ItP 

0 

(7V-1) 

4(iV + l) 
N-3 

4 {N + 2) 
N-3 

4 (AT + 2) 
N-3 

2N 

0 

4{N -1) 

4{N -2) 

4iN-2) 

2N 


Table 2: Matrix C 


Y, (-1)^^' A, D, C] = (-1)^* (- A? + A./2) A (C.26) 

j 

were we introduced the symmetry index Sj = s[7^j] via 


piab _ ^_ p^Si piba _ ^_ p^Si piab 


(C.27) 


For symmetric representations Sj = 0 while for antisymmetric ones = 1 (table |^). The matrix C 
arises when a projector multiplies a product of two matrices T: 


pi fSd 

a7 


rpa rpa 

aa d5 


= N Xj = N i-ir pV ^ P'aS 


jjad _ 


ad 


— JSJ-(^_iyi+Sk pk05 


Q7 


(C.28) 


We also have 

P%' P'li Tfid Tib = ^ A, , (C.29) 

We now list some useful properties obeyed by the tensors Z. Under the exchange of indices they 
behave as: 


Z+ 


ab 

cd 


fy-\- cd , 
^ ab ’ 



y— cd . 
^ ab ’ 


/ 7 + ab 
^ cd 


ry— ab . 

^ dc ’ 


ry-\- ab 

^ cd 


(C.30) 


Multiplication properties of Z are 


ab 
^ cd 

PlSAlii : 

= ±i(Z+ - Z-)ti-, 





ab 
^ cd 


= i(Z+ + Z-)ti 





/vi ab 
^ cd 

y^ cd 
^ kl — 

T (IV^ - 4) (P|8.4] + 





ry-\- ab 
^ cd 

^ kl — 

(iV2 _ 4) (p[8^] _ p[8g])-^ 




r/-\- d5 
^ aa 

yd yd 

^a'y ^0d - 

= z-llNxy^;^^ = 


pioU“ \ 

r )da 7 

pjad _ 
^ 'r0 - 

0 

ry— dS 

^ aa 

ya yd 

^a'y ^ 0d 

= Z^llNX,P0^ = 


\ pj ad _ 

7^7/3 - 

: -NZ 

+ (Sa \ 
7/3 -^8 

Z±0S 

a'y 

T^i bd yd 
^ aa ^60 

1 

0 

0 

II 





Z±05 

a7 

pi d0 pa 
^ ca ^ a 7 

Tl, = 

- 





2 


(C.31) 
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